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1  Executive  Summary 


This  project  was  initiated  as  a  response  to  an  RFP  soliciting  novel  and  innovative  ap¬ 
proaches  to  optical  power  limiting.  The  central  thrust  of  the  original  contract  was  the 
investigation  of  liquid  crystal  systems,  including  pure  and  composites,  as  potential  sources 
for  novel  nonlinear  optical  effects  that  could  be  utilised  in  optical  power  limiting  applica¬ 
tions. 

The  initial  phase  of  our  activities  investigated  polymer  dispersed  liquid  crystal  systems 
in  active  and  passive  OPL  configurations.  Polymer  dispersed  liquid  crystals,  (PDLCs) 
had  previously  been  shown  to  be  effective  optical  shutters  with  essentially  no  insertion 
loss  due  to  absorption  and  thus  “water  clear”  optical  properties  in  the  transmitting  state. 
PDLCs  can  be  switched,  through  electro-optic  mechanisms,  to  essentially  opaque  (milky 
white)  states  in  the  scattering  mode.  We  assessed  our  findings  as  exceptionally  effective 
for  OPL  with  the  active  electro-optic  configuration  and  promising  for  passive  operation. 

The  later  effort  of  the  contract  followed  the  directions  of  interest  indicated  by  DARPA 
and  the  COTR,  seeking  results  for  OPL  using  liquid  crystals  in  the  nanosecond  regime. 
A  smaller  effort  continued  in  the  CW  and  picosecond  regime.  With  interrupted  funding 
and  a  reduced  rate  of  funding,  the  effort  continued  from  June  1988  until  June  1993. 

A  broad  selections  of  liquid  crystals  were  synthesized  or  purchased  as  part  of  a  screening 
program  in  nonlinear  materials  properties.  The  nonlinear  optical  properties  were  mea¬ 
sured  using  the  Z-scan  technique  to  determine  the  nonlinear  indices  of  refraction  and  the 
nonlinear  absorption  coefficients.  These  properties  were  determined  for  CW,  nanosecond 
and  picosecond  lasers.  Materials  were  also  tested  for  optical  power  limiting  performance. 
Several  good  candidates  were  found  with  large  nonlinearities  and  good  limiting  charac¬ 
teristics. 


Our  group  was  the  first  to  discover  5CB  as  a  nonlinear  liquid  crystal  with  OPL  po¬ 
tential.  5CB  (4-cyano-4’-n-pentylbiphenyl),  8CB  (4-cyano-4’-n-octylbiphenyl)  and  a  few 
selected  few  compounds  were  investigated  in  detail,  while  approximately  30  other  com¬ 
pounds  were  evaluated.  We  carried  out  an  effort  to  identify  the  nonlinear  mechanisms  as 
well  as  screen  promising  materials.  The  presence  of  two-photon  absorption  and  fifth  order 
nonlinearities  were  determined  present  as  important  mechanisms  in  the  best  materials. 


In  addition  to  the  nonlinear  studies,  a  separate  subproject  involving  Unear  properties 
of  polymer  dispersed  liquid  crystals  in  the  infrared  region  was  performed.  The  modula¬ 
tion  properties  of  these  materials  were  found  promising,  but  difficulties  were  identified 
involving  insertion  loss  and  response  time.  ,  By. 
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2  Introduction 


This  research  effort  carried  out  a  wide  range  of  investigations  of  both  linear  and  nonlinear 
optical  phenomena  in  liquid  crystals  in  order  to  evaluate  their  effectiveness  in  certain 
optical  power  limiting  and  optical  modulation  applications.  The  areas  of  investigation 
included: 


1.  CW  Ar+  power  limiting  by  liquid  crystal  composites 

2.  Nanosecond  Nd  YAG  laser  pulse  limiting  by  pure  liquid  crystal 

3.  Picosecond  Nd  YAG  laser  pulse  limiting  by  pure  liquid  crystal 

4.  IR  shuttering  capabilities  of  polymer  liquid  crystal  composites. 

These  activities  were  supported  by  efforts  in  three  areas: 


•  Experimental  laser  measurements 

•  Materials  preparation  and  synthesis 

•  Theoretical  analysis. 


The  body  of  this  report  outlines  the  specific  research  efforts  and  summarises  the  results. 
Additional  details  may  be  found  in  the  Quarterly  Scientific  and  Technical  Reports  which 
were  filed  with  the  COTR  during  the  course  of  investigation. 
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3  CW  Regime 


The  main  thrust  of  the  first  stage  of  the  experimental  effort  was  the  investigation  of  the 
possibilities  for  power  limiting  of  the  polymer  dispersed  Hquid  crystal  (PDLC)  films  with 
continuous  wave  (CW)  lasers.  This  proceeded  in  accordance  with  the  proposed  work 
by  investigating  both  passive  and  active  electronic  feedback  configurations.  The  effort 
resulted  in  an  active  feedback  system  meeting  the  requirements  for  an  Optical  Power 
Limiting  (OPL)  device  in  the  CW  regime,  i.e.  transparent  at  low  intensity  and  blocking 
high  intensity  transmission. 


3.1  Active  Linear  Feedback  Optical  Power  Limiting  for  CW 
Ar+  laser 

The  active  feedback  investigations  studied  a  combined  device  design  with  an  electronic 
sensor  that  controlled  a  light  shutter  using  liquid  crystals.  The  light  shutter  technology 
is  based  on  the  polymer  dispersed  liquid  crystal  (PDLC)  technology  established  in  liquid 
crystal  display  (LCD)  applications. 

The  use  of  PDLCs  in  display  applications  is  a  mature  technology.  It  is  based  on  scat¬ 
tering  of  light  by  liquid  crystal  inclusions  in  a  transparent  polymer.  While  the  materials 
and  geometry  of  a  PDLC  may  vary  greatly,  there  is  a  common  banc  structure  from  which 
all  applications  derive.  This  basic  structure  is  a  dear  polymer  which  has  impregnated  (by 
a  variety  of  processes)  spherical  inclusions  of  liquid  crystals.  (Figure  3  shows  this  geome¬ 
try  for  example.)  The  liquid  crystal  inclusions  are  essentially  nonabsorbing,  so  when  the 
index  of  refraction  of  the  liquid  crystal  is  the  same  as  the  polymer,  then  the  i, 

“water  clear”.  Because  all  liquid  crystal  materials  are  birefringent,  they  have  two  differ¬ 
ent  indices  of  refraction  depending  on  the  orientation  of  the  liquid  crystal  relative  to  the 
direction  and  polarization  of  light.  This  direction  of  orientation  is  electrically  controlled. 
This  control  forms  the  basis  of  the  active  feedback  mechanism  investigated. 

The  PDLC  films  studied  consisted  of  the  hquid  crystal  E7  dispersed  in  the  epoxy  matrix 
made  using  Epon  282,  MK107  and  Capcure.  The  films  can  be  switched  from  a  scattering 
opaque  state  to  a  dear  transparent  state  by  the  application  of  a  voltage  across  the  film. 
We  constructed  a  feedback  circuit  so  that  the  voltage  applied  across  the  film  was  a  function 
of  the  intensity  of  light  transmitted  by  the  film.  The  optical  power  transmitted  by  the 
sample  was  measured  by  a  photodiode;  the  output  voltage  of  the  photodiode  amplifier 
circuit  was  used  as  the  input  of  a  linear  feedback  circuit.  The  amplitude  of  the  1kHz 
sinusoidal  voltage  applied  to  the  film  was  therefore  a  linear  function  of  the  transmitted 
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Figure  1:  Schematic  of  experimental  arrangement  for  active  feedback  OPL. 
power.  The  schematic  diagram  of  the  experimental  arrangement  is  shown  in  Fig.  1. 

The  transmittance  of  the  film  depends  on  the  applied  voltage,  and,  in  the  presence  of 
feedback,  becomes  a  function  of  incident  power.  The  response  of  the  sample  in  the  case  of 
negative  feedback  exhibits  optical  power  limiting  behavior  as  shown  in  Fig.  2  for  different 
feedback  parameters. 

In  the  case  of  negative  feedback,  for  some  range  of  optical  input  power,  the  system  has 
two  stable  and  one  unstable  states  for  a  given  input.  This  optical  bistability  may  be  of 
interest  for  optical  memory  applications. 

The  OPL  characteristics  of  this  feedback  arrangement  are  excellent  and  offer  promise 
for  use  in  CW  protection  applications.  There  is  considerable  latitude  in  the  selection 
of  electric  circuit  parameters,  such  as  driving  voltage  and  sensor  sensitivity.  Similarily,  a 
variety  of  liquid  crystal  and  polymer  materials  may  be  used  and  the  liquid  crystal  droplets 
can  be  varied  in  size  and  concentration.  Because  of  this  flexibility,  the  system  parameters 
may  be  optimized  for  a  variety  of  applications. 

The  work  on  this  aspect  of  the  project  was  designed  to  provide  a  demonstration  of  a 
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Figure  2:  Optical  power  limiting  response  for  different  feedback  parameters. 

CW  OPL  device  concept.  A  prototype  of  this  device  was  delivered  to  the  COTR.  Work 
was  completed  on  this  phase  of  the  project.  Farther  efforts  could  be  directed  toward 
meeting  specific  device  requirements  if  such  requirements  are  provided.  In  particular, 
response  times  and  dynamic  range  could  be  tested  and  optimized  in  the  future. 


3.2  Passive  OPL  Study  of  PDLCs 

One  of  the  contract  requirements  was  to  investigate  passive  OPL  properties  of  liquid 
crystal  over  a  wide  range  of  time  scales.  Passive  OPL  has  inherent  advantages,  such  as 
not  relying  on  electronic  circuits,  but  developing  materials  properties  and  geometries  is 
inherently  more  difficult.  For  CW  laser  beams,  active  feedback  devices  such  as  that  de¬ 
scribed  in  the  previous  section  may  be  expected  to  be  superior.  Nonetheless,  passive  OPL 
behavior  was  studied  and  characterized  in  PDLCs  as  required.  Two  sets  of  experiments 
were  carried  out  to  study  passive  optical  switching  in  PDLCs.  The  mechanisms  available 
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for  CW  twitching  included  heating  effects  and  optical  field  induced  reorientation  of  the 
liquid  crystal  molecules  in  the  droplets. 

The  first  experiments  studied  these  two  switching  mechanisms  and  the  geometry  used 
in  the  study  was  an  initially  transparent  PDLC  film.  The  observation  of  switching  is  then 
observed  as  induced  scattering  by  the  liquid  crystal  droplets  as  illustrated  schematically 
in  Fig.  3. 

The  geometry  for  optical  switching  starts  from  the  “on”  or  transparent  state  described 
above  with  a  quasi-static  field  applied  perpendicular  to  the  plane  of  the  film.  The  optical 
field  of  an  argon-ion  laser  is  then  incident  normally  on  the  film  and  the  electric  field  of 
the  light  is  perpendicular  to  the  quasi-static  field  maintaining  the  transparent  “on”  orien¬ 
tation  of  the  droplets.  In  theory,  when  the  electric  field  of  the  light  becomes  comparable 
to  the  quasi-static  field,  the  competition  will  cause  the  loss  of  orientation  or  possibly 
reorientation  of  the  droplets.  We  found  that  the  intensity  of  our  10  watt  Argon  ion  laser 
was  adequate  to  this  purpose. 


Figure  3:  Light  scattering  results  when  the  optical  field  induces  a  change  in  the  dielectric 
properties  of  the  inclusions  producing  optical  index  mismatch. 
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We  set  up  this  experiment  with  a  chopped  Argon  ion  beam  incident  on  an  initially 
transparent  PDLC.  The  transmittance  of  the  PDLC  was  monitored  by  a  small  He-Neon 
laser.  The  layout  of  the  apparatus  is  shown  in  Fig.  4. 

At  low  intensity  of  the  Argon  ion  laser,  the  probe  beam  is  essentially  nnattennated  (90% 
transmittance).  At  high  intensity,  we  found  that  transmission  was  strongly  attenuated. 
Fig.  5  shows  transmittance  measured  by  the  probe  beam  as  function  of  the  At  pump 
beam  intensity. 

The  maximum  decrease  in  probe  beam  intensity  is  about  80%  and  is  dependent  on  the 
relative  orientation  of  polarisation  of  the  pump  and  probe  beams.  This  difference  in  the 
polarisation  indicates  that  both  optical  reorientation  and  thermal  effects  are  present. 

A  second  set  of  experiments  on  laser  inducing  scattering  in  PDLCs  was  performed  to 
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farther  elucidate  the  competition  between  the  optical  field  reorientation  and  the  heating 
effect  in  the  switching  experiments  on  PDLCs.  There  are  two  essential  differences  between 
the  data  taken  in  the  second  set  of  experiments.  Firstly,  PDLCs  were  artificially  heated 
to  temperatures  approaching  the  nematic-isotropic  transition.  Secondly,  we  performed 
experiments  on  samples  with  no  stabilizing  static  field  applied  to  the  droplets  initially. 
We  were  therefore  trying  to  switch  from  a  scattering  state  to  a  transmitting  state  in  these 
experiments. 


A r+  Loser  Intensity  [kW/cm^J 


Figure  5:  Minimum  transmittance  as  a  function  of  pump  beam  intensity. 
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Transmittance  Increase  (%) 


A>  expected,  there  is  no  effect  whatever  at  room  temperature.  At  higher  temperatures, 
the  elastic  and  viscous  forces  in  the  liquid  crystal  droplets  weaken  and  we  expect  an 
optical  held  to  be  more  effective  in  reorienting  liquid  crystal  molecules  into  a  transmitting 
geometry.  This  is  indeed  what  was  observed  as  is  shown  in  Fig.  6. 

The  particular  importance  of  these  results  lies  not  in  the  use  of  the  effect  (although 
it  is  not  small)  but  rather  in  the  demonstrated  optical  orientability  of  a  random  PDLC. 
To  our  knowledge,  this  was  the  first  such  demonstration  although  we  have  had  informal 
reports  of  attempts  by  researchers  at  Hughes,  University  of  California  at  Berkeley [1],  and 
Celanese  in  similar  systems. 
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Figure  6:  Optical  field  induced  transmission. 
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3.3  Nonlinear  Materials  Parameters  in  the  CW  regime 

Pate  liquid  crystal  materials  were  studied  extensively  is  addition  to  the  PDLC  systems. 
Results  tabulating  CW  and  pulsed  materials  characteristics  are  contained  in  Appendix  A 
and  Sec.  5.1.  Here  we  report  the  experiments  on  selected  compounds. 

We  carried  out  systematic  measurements  of  the  nonlinear  refractive  indices  n*  of  pure 
liquid  crystal  materials.  The  method  of  measurement  is  the  c-scan  technique,  developed 
by  E.  van  Stryland  et  al.  at  CREOL.  This  is  essentially  a  measurement  of  self-focusing  by 
the  sample;  the  nonlinear  refractive  index  of  the  sample  can  be  readily  determined  from 
intensity  of  the  transmitted  radiation  in  the  far  field.  Since  liquid  crystals  are  birefringent, 
we  can  measure  two  nonlinear  refractive  indices  (and  the  nonlinear  birefringence). 

During  the  z-scan  measurements,  the  samples  are  translated  along  the  direction  of 
beam  propagation  across  the  focal  point  of  a  beam  focusing  lens.  The  samples  are  ho¬ 
mogeneously  aligned  (parallel)  monodomain  samples  of  the  liquid  crystal  5CB,  between 
parallel  glass  plates.  The  sample  thickness  varied  from  25  pm  to  120pm.  The  Ar+  beam 
was  chopped  by  an  electromechanical  shutter  to  produce  pulses  with  a  width  of  approxi¬ 
mately  1G  ms.  The  far-field  transmittance  was  measured  with  the  laser  polarization  both 
parallel  and  perpendicular  to  the  nematic  director.  (The  nematic  director  is  the  direction 
the  molecules  point,  and  hence  the  direction  along  which  the  polarization  of  the  light 
experiences  the  greatest  linear  index  of  refraction.) 

Typical  results  for  the  z-scan  measurements  of  transmission  intensity  versus  position  of 
the  sample  are  shown  in  figures  7  and  8.  The  figure  shows  the  results  of  two  measurements, 
one  with  the  polarization  parallel  to  the  director  (Fig.  8)  and  one  with  the  polarization 
perpendicular  to  the  director  (Fig.  7).  The  sample,  5CB,  has  a  very  typical  structure  of 
a  liquid  crystal  (see  structures  in  Appendix  A). 

A  sophisticated  fitting  procedure  was  employed  to  determine  the  nonlinear  optical 
parameters  of  the  sample  material  from  the  z-scan  measurements.  Even  without  such 
analysis,  it  can  be  seen  from  the  above  figure  that  the  nonlinear  response  of  the  liquid 
crystal  is  extremely  large.  First,  observe  that  the  sign  of  the  nonlinear  index  of  refrac¬ 
tion  is  opposite  for  the  two  polarization.  As  the  sample  is  translated  the  perpendicular 
polarization  intensity  first  decreases  then  increases.  The  opposite  is  true  for  the  parallel 
polarization.  Secondly,  notice  that  the  decrease  in  intensity  reaches  a  minimum  of  around 
10%  of  the  unaltered  intensity.  This  sample  is  exhibiting  optical  power  limiting  and  it  is 
quite  a  large  effect.  The  effect  is  similar  in  other  liquid  crystals  with  similar  molecular 
structures. 

The  quantitative  analysis  of  our  results  with  the  Ar+  laser,  at  A  =  514nm,  with  10ms 
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Figure  7:  (left).  Z-scan  of  25  pm  thickness  sample  of  5CB  with  polarisation  perpendicular 
to  director,  CW  AR+  laser  at  514  nm. 

Figure  8:  (right).  Z-scan  of  25  pm  thickness  sample  of  5CB  with  polarisation  parallel  to 
director  using  CW  AR+  laser  at  514  nm. 
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pulses,  yielded  the  following  values  for  the  nonlinear  index  of  refraction  of  5GB: 


•  for  polarisation  parallel  to  the  director  nj  w  —  10“a  esu, 

•  for  polarisation  perpendicular  to  the  director  nj  »  10~2  esu. 


We  see  then  that  the  nonlinearity  is  large,  bnt  of  the  opposite  sign  than  is  expected 
from  optical  field  reorientation.  In  particular,  the  negative  nj  means  the  material  is  *e(f 
defocuaing  for  polarisation  parallel  to  the  director.  This  is  suggestive  of  large  thermal 
effects  which  were  also  investigated  and  are  described  in  Sec.  3.3.1. 


3.3.1  Thermal  Behavior  in  the  CW  Regime 


The  negative  value  of  the  nonlinear  index  in  the  direction  of  molecular  ordering  suggests 
that  heating  may  be  the  operative  mechanism.  In  such  an  event,  the  increase  in  tem¬ 
perature  would  decrease  the  order  parameter,  and  the  parallel  index  of  refraction  would 
decrease  as  well.  The  opposite  would  be  true  for  the  index  of  refraction  when  the  polar¬ 
isation  is  perpendicular  to  the  direction  of  orientation. 

The  physical  properties  of  liquid  crystals  change  dramatically  as  they  are  heated  toward 
the  isotropic  (normal  liquid)  phase.  The  above  results  are  for  a  temperature  a  few  degrees 
below  the  melting  temperature.  To  determine  the  change  of  the  nonlinearity  as  a  function 
of  temperature,  the  temperature  dependence  of  these  quantities  was  also  determined.  A 
divergence  of  the  nonlinear  index  as  the  isotropic  transition  approached  is  shown  in 
Fig.  9. 

The  temperature  dependence  in  Fig.  9  is  for  the  nematic  liquid  crystal  5CB  using 
10ms  Ar+  laser  pulses  at  A  =514nm.  The  nonlinear  birefringence  diverges  as  temperature 
approaches  the  nematic-isotropic  transition  temperature.  The  extremely  small  nonlinear¬ 
ity  above  the  transition  temperature  shows  that  it  is  the  properties  of  the  liquid  crystal 
nematic  phase  which  is  the  origin  of  the  nonlinearity  and  not  a  property  inherent  in  the 
individual  molecules.  More  accurately,  the  molecules  ate  not  able  to  contribute  to  the 
optical  nonlinearity  effectively  unless  they  are  in  the  liquid  crystal  phase. 
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This  it  compelling  evidence  for  a  strong  thermal  component  in  the  nonlinear  response 
in  the  CW  regime.  Nonetheless,  there  it  alto  dear  evidence  for  optical  reorientation  in 
the  PDLC  results  of  Sec  3.2.  (See  Sec.  7.2  for  a  related  theoretical  ditcutskm.) 


Figure  9:  The  temperature  dependence  of  ftyj  and  n2x 
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4  Nanosecond  Nd  YAG  Nonlinear  Optical  Response 


A  systematic  study  of  nonlinear  parameters  of  a  variety  of  liquid  crystals  was  carried  out 
using  our  Q-switched  Nd:YAG  laser,  with  nanosecond  pulses.  Typical  results  for  5CB  for 
polarisations  parallel  and  perpendicular  to  the  director  are  shown  in  Figs.  10  and  11.  The 
effect  for  the  parallel  polarisation  is  large,  for  the  perpendicular  polarisation  is  quite  small. 
The  third-order  hyperpolarisability  in  the  nanosecond  regime  is  «  10~9  esu,  some  two 
orders  of  magnitude  greater  than  that  of  CS2.  Results  for  8CB,  which  exhibits  a  smectic 
phase,  are  considerably  greater.  Theoretical  work  provided  models  for  the  mechanisms 
for  these  effects  and  is  discussed  in  Sec.  7.  Additional  work  on  measurements  of  other 
liquid  crystalline  materials,  and  on  molecular  structure/physical  property  relationships  is 
tabulated  in  Appendix  A  and  in  Sec.  5.1  where  we  discuss  the  mechanisms  involved  in 
the  nanosecond  nonlinearities. 

The  measurements  on  the  5CB  were  used  to  determine  the  nonlinear  index  of  refraction 
and  nonlinear  absorption  coefficient  at  A  =532nm.  We  obtained  the  values: 

•  n2  =  —54  x  10~n  esu  for  polarization  parallel  to  director. 

•  n2  =  +8.3  x  10~u  esu  for  polarization  perpendicular  to  director. 

•  »j  =  -24  x  10"n  esu  for  the  isotropic  phase  . 


The  nonlinear  absorption  coefficient  0  has  also  been  measured  with  the  results: 


•  0  ~  265  cm/GW  for  polarization  parallel  to  director. 

•  0  =  36  cm/GW  for  polarization  perpendicular  to  director. 

•  0  =  114  cm/GW  for  for  the  isotropic  phase. 

The  key  result  is  that  n2  for  5GB  in  the  nematic  phase  in  the  nanosecond  regime  for 
the  wavelength  A  =  532nm  is  nearly  two  orders  of  magnitude  greater  than  CS2.  We  have 
verified  the  accuracy  of  our  measurements  by  measuring  n2  for  CS2,  and  have  obtained 
values  in  agreement  with  those  in  the  literature. 

We  have  also  attempted  to  measure  the  nonlinear  response  at  the  fundamental  wave¬ 
length  A  =  1.06/im  of  the  Q-switched  Nd:YAG  laser.  Our  experimental  results  for  5CB 
show  very  weak  nonlinearities  at  this  wavelength,  up  to  the  damage  threshold  of  the  sam¬ 
ple  cell  (not  the  liquid  crystal).  (For  a  discussion  of  damage,  see  Sec.  6.)  The  nonlinear 
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Figure  10:  (left).  Z-scan  of  120  pm  thickness  sample  of  5CB  with  polarisation  parallel  to 
director  using  8ns  pulses  at  532  nm. 

Figure  11:  (right).  Z-scan  of  120  pm  thickness  sample  of  5CB  with  polarisation  perpen¬ 
dicular  to  director  using  8ns  pulses  at  532  nm. 

contribution  to  the  transmittance  could  not  be  effectively  separated  from  noise  because 
of  its  small  magnitude. 

As  in  the  CW  case,  we  see  then  that  the  nonlinearity  is  large,  but  of  the  opposite 
sign  than  is  expected  from  optical  field  reorientation.  In  particular,  the  negative  n2 
means  the  material  is  self  defocvsing  for  polarisation  parallel  to  the  director.  In  the  CW 
measurements,  we  interpreted  the  negative  value  of  n2  for  parallel  polarisations  to  be 
due  to  heating.  In  the  nanosecond  regime  this  would  be  a  very  peculiar  result  because 
conventional  wisdom  would  suggest  that  thermal  effects  are  too  slow  to  be  operative  in  the 
nanosecond  regime.  The  temporal  characteristics  of  the  nonlinear  response  were  therefore 
experimentally  investigated  and  are  reported  in  Sec.  4.1.  (The  possible  role  of  thermal 
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effects  were  also  investigated  theoretically  and  are  discussed  in  Sec.  7.2.) 


4.1  Time  Dependence  of  Nonlinear  Response:  Two-Pulse  Mea¬ 
surements 

In  order  to  identify  the  dominant  mechanism  responsible  for  the  observed  large  nonHn- 
earities  in  nanosecond  regime,  it  is  useful  to  have  an  estimate  of  the  characteristic  time  of 
the  contributing  process.  To  estimate  this  time,  we  have  devised  a  time  resolved  double 
pulse  Z-scan  scheme.  The  experiment  setup  is  shown  in  Fig.  12. 


Figure  12:  Experimental  set  up  for  two-pulse  measurements 
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Figure  13:  Transmitted  pulses  when  the  sample  was  placed  far  from  the  focal  point 

In  this  experiment,  a  portion  of  the  beam  from  the  laser  is  first  delayed  and  then 
recombined  with  the  original  beam,  resulting  in  two  pulses  of  nearly  identical  amplitude 
separated  in  time.  By  performing  Z-ican  measurements  using  these  two  pulses,  it  is 
possible  to  determine  if  the  first  pulse  affects  the  response  to  the  second.  If  it  does,  then 
the  decay  time  of  the  process  can  not  be  less  than  the  temporal  separation  of  the  pulses. 
For  this  scheme  to  work,  it  is  essential  that  at  the  sample  the  two  beams  be  parallel  and 
that  both  pulses  pass  through  the  same  region  of  the  sample,  a  nontrivial  task. 

\ 

The  time  delay  between  the  two  pulses  can  be  varied  by  changing  the  length  of  the 
delay  line.  The  energy  of  each  pulse  can  be  a4justed  independently  by  changing  the  angles 
between  the  half  wave  plate  and  the  Gian  prisms.  In  our  experiment,  the  two  pulses  had 
nearly  the  same  energy  and  a  time  delay  down  to  llns.  The  incident  and  transmitted 
intensities  were  recorded  on  the  oscilloscope.  The  sample  was  5CB  at  24  degrees  C.  The 
incident  polarization  of  both  pulses  was  parallel  to  the  director.  Figs.  13  and  14  show  the 
oscilloscope  record  when  a  small  aperture  was  placed  in  front  of  the  photodetector  PD1. 
In  this  geometry,  a  nonlinear  response  removes  energy  from  the  beam. 

Fig.  13  shows  the  transmitted  pulses  when  the  sample  was  placed  far  from  the  focal 
point.  In  this  geometry,  the  intensity  in  the  sample  is  low  and  hence  no  nonlinear  response 
is  expected.  As  expected,  Fig.  13  shows  that  the  two  transmitted  pulses  have  almost  the 
same  height.  This  confirms  that  the  intensity  is  low  and  nonlinear  effects  are  negligible. 
However,  when  the  sample  is  moved  dose  to  the  focal  point,  at  a  particular  sample 
position,  the  first  pulse  is  suppressed  due  to  nonlinear  absorption  and  nonlinear  refraction, 
while  the  second  one  is  increased  rather  than  decreased.  Such  an  example  is  shown  in 
Fig.  14. 
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Figure  14:  Transmitted  pulses  when  the  sample  was  placed  near  the  focal  point 

The  increase  in  the  amplitude  of  the  second  transmitted  pulse  is  clearly  due  to  the 
effect  of  the  first  pulse,  since  when  the  first  pulse  is  blocked,  the  second  pulse  decreases 
also.  This  behavior  implies  that,  regarding  nonlinear  refraction,  the  second  pulse  “sees” 
the  effect  of  the  first  pulse.  However,  if  the  total  transmittance  is  measured  without  an 
aperture  in  front  of  the  detector,  i.e.  all  transmitted  light  is  detected,  the  two  transmitted 
pulses  have  the  same  response  regardless  of  the  sample  position  relative  to  the  focal  point. 
This  implies  that,  regarding  nonlinear  absorption,  the  second  pulse  does  not  see  the  effect 
of  the  first. 

These  observations  raised  significant  questions  about  the  role  of  orientational  order 
in  nonlinear  response.  If  the  first  pulse  changed  the  orientational  character  of  medium, 
and  the  nonlinear  absorption  were  anisotropic,  then  the  absorption  should  change  for  the 
second  pulse  as  well.  One  of  the  ways  to  change  orientational  order  in  liquid  crystals 
is  to  change  the  temperature,  since  they  become  more  ordered  at  low  temperature.  Ex¬ 
periments  on  the  temperature  dependence  of  nanosecond  nonlinearities  were  accordingly 
performed  and  are  described  in  Sec.  4.2. 


4.2  Temperature  Dependence  of  Nanosecond  Nonlinear  Re¬ 
sponse 

We  have  also  carried  out  measurements  of  the  temperature  dependence  of  the  nonlinear 
response  of  5CB  on  the  nanosecond  scale  at  532nm  wavelength.  We  found,  as  in  millisec* 


ond  case,  that  najj  and  n,x  in  the  nematic  phase  hare  opposite  signs,  and  that  there  » 
an  abrupt  change  in  these  values  across  the  nematic-isotropic  transition.  However,  unlike 
in  the  millisecond  case,  the  values  of  »aj|  and  njj.  do  not  change  appreciably  with  tem¬ 
perature  in  the  nematic  phase,  and  remain  nearly  constant  up  to  the  nematic-isotropic 
transition  temperature. 

Figs.  15  undid  show  the  temperature  dependence  of  the  parallel  component  of  and 
0.  We  note,  as  shown  in  Fig.  15,  that  njj)  is  constant  to  within  5%  from  0  to  10  degrees 
C  below  the  transition;  the  transition  was  approached  to  within  less  than  400mK.  Similar 
behavior  holds  for  0\\  in  Fig.  16.  However,  there  is  an  abrupt  jump  at  the  transition 
temperature  both  for  n2\\  and  0\\. 


Figure  15:  Temperature  dependence  of  the  parallel  component  of  n2. 
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Figure  16:  Temperature  dependence  of  the  parallel  component  of  0 
4.3  Higher  Order  Nonlinearities 

We  have  also  measured  the  intensity  dependence  of  the  nonlinear  absorption  coefficient  /fy 
and  refractive  index  n2||  for  5CB  at  532nm  using  ns  pulses.  The  sample  has  a  thickness  of 
25 pm  and  the  sample  temperature  is  24  degrees  C.  In  this  set  of  measurements,  for  each 
incident  energy,  we  performed  a  Z-scan  measurement  and  obtained  the  nonlinear  refraction 
and  absorption  coefficients.  In  this  way,  n2  and  were  determined  as  function  of  incident 
intensity.  Fig.  18  shows  the  result  for  the  geometry  where  the  incident  polarisation  is 
parallel  to  the  director.  The  horisontal  axis  represents  on-axis  laser  intensity  I0,  defined 
as  Jo  =  where  P  is  the  power  and  Mo  is  the  beam  waist. 

As  can  be  seen  in  Fig.  17,  the  nonlinear  refractive  index  tty)  remains  essentially  con¬ 
stant  within  our  experimental  error,  as  expected  for  a  third-order  process.  Interestingly, 
however,  the  effective  nonlinear  absorption  increases  linearly  with  the  intensity.  This 
has  implications  for  the  mechanisms  operating  to  produce  the  nonlinearity.  We  discuss 
these  mechanisms  further  in  Secs.  5.1  and  7.3,  but  we  note  that  the  increase  in  nonlinear 
absorption  or  refractive  index  with  intensity  serves  to  enhance  the  OPL  behavior. 
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Figure  17:  Intensity  dependence  of  the  nonlinear  refractive  index  n2||  for  5CB  at  532nm 
using  7  ns  pulses 
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Figure  18:  Intensity  dependence  of  the  nonlinear  absorption  coefficient  /fy  for  5CB  at 
532nm  using  7  ns  pulses 
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4.4  Photoacoustic  Measurements 


We  also  attempted  photoacoustic  measurements  where  an  attempt  was  made  to  detect 
an  acoustic  wave  generated  by  a  nonlinear  pulse  under  conditions  of  strong  nonlinear 
absorption.  We  did  not  have  an  unambiguous  answer  to  the  question  of  where  the  energy 
went  which  was  removed  from  the  beam  during  OPL  circumstances,  thus  we  sought 
to  investigate  sound  generation  as  a  process  for  dissipation  of  the  absorbed  energy.  A 
sensitive  accelerometer,  consisting  of  a  piezoelectric  element,  was  attached  to  one  glass 
plate  of  a  120  /im  thick  SCB  sample.  The  sample  was  then  exposed  to  532nm  pulses.  The 
incident  polarization  was  parallel  to  the  director  and  the  pulse  energy  was  64  pJ.  With  the 
sample  positioned  near  the  focal  point  of  the  lens,  the  transmittance  was  about  0.5.  At  this 
intensity,  a  signal  was  detected  by  the  accelerometer,  indicating  that  acoustic  excitations 
play  do  play  some  role  in  the  dissipation  process.  Unfortunately,  during  the  duration  of 
this  effort,  we  were  unable  to  refine  these  measurements  to  obtain  a  quantitative  measure 
of  the  amount  of  energy  dissipated. 


Figure  19:  Optical  power  limiting  behavior  of  25  pm  sample  of  5CB. 

4.5  Nanosecond  Optical  Power  Limiting 

Optical  power  limiting  behavior  in  these  materials  on  the  nanosecond  time  regime  has 
also  been  observed.  In  Fig.  19  we  show  the  OPL  response  of  a  25  pm  thick  sample  of 
nematic  5CB  (homogeneous  alignment)  to  8ns  pulses  at  532nm. 

Notice  that  the  shape  of  the  curve  in  Fig.  19  shows  a  nearly  ideal  OPL  behavior.  The 
energy  threshold  here  is  «  lmJ  with  a  500pm  spot  sise.  This  threshold  is  greatly  reduced 
in  geometries  appropriate  to  devices  used  in  the  field.  (See  Sec.  4.5.1.)  With  a  faster 
lens  to  give  a  10pm  spot  sise,  this  threshold  would  be  reduced  to  0.4pJ.  We  produce  this 
particular  plot  here  because  believe  this  was  the  first  demonstration  among  any  of  the 
contractors  of  effective  OPL  by  nematic  liquid  crystals  in  the  nanosecond  regime. 


4.5.1  ARDEO  DVO  Test  Bed 

Under  this  contract  we  provided  two  "best  samples’*  for  evaluation  with  the  ARDEO 
DVO  test  bed.  These  samples  were  of  8CB  and  differed  in  thickness  and  “doping”,  with 
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dye  (anthraquinone  dye  dm  from  BDH  it  itself  a  liquid  crystal).  Sample  1  is  240  mm  thick 
and  has  a  dye  concentration  of  0.05%.  Sample  2  is  120  mm  thick  with  dye  concentration 
0.02%.  These  are  very  thick  samples  for  optimum  effect  of  the  nonlinear  response.  The 
nominal  transmission  of  our  “best  sample0  was  80%,  in  keeping  with  the  transmission 
standards  required. 

Prior  to  delivery,  we  replicated  the  required  optics  and  performed  our  own  measure¬ 
ments  of  OPL.  These  samples  were  also  characterised  for  nonlinear  absorption  and  non¬ 
linear  refraction. 

Our  measurements  were  performed  with  a  “top  hat”  beam  profile  as  specified  for  use  in 
the  test  bed.  Fig.  20  shows  the  beam  profile.  It  is  “dirty”  as  required,  showing  significant 
spikes  across  the  profile.  This  is  a  characteristic  expected  under  field  conditions. 


2.5  5.0 

position  law] 


Figure  20:  Beam  profile  used  in  measuring  samples  provided  for  test  bed 
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Figure  21:  Z-scan  to  determine  nonlinear  refractive  index  using  top  hat  beam 

Samples  were  characterized  ««ng  the  Z-scan  technique  for  nonlinear  index  and  nonlin¬ 
ear  absorption  as  shown  in  Figs.  21  and  22  respectively. 

Doped  samples  of  8CB  give  the  largest  nonlinear  index  of  refraction  and  the  largest 
nonlinear  index  of  refraction  of  any  materials  we  have  tested.  As  stated  above,  80% 
transmission  was  maintained  with  a  doping  of  leas  than  0.05%  by  weight.  The  values 
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Figure  22:  Z-scan  to  determine  nonlinear  absorption  using  top  hat  beam 
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n3jj  =  —152  x  10~neau  and  njx  —  +10.5  x  10~nejn  were  obtained  for  polarisations 
parallel  and  perpendicular  to  the  director.  In  the  isotropic  phase,  the  nonlinear  index, 
n2l  =  —48  x  10-n  eiu  was  obtained.  Compared  to  the  pure  material,  the  doping  increased 
the  largest  nonlinearity  by  a  factor  of  six. 

Large  values  of  the  nonlinear  absorption  were  also  obtained.  The  values  0\\  =  380 
cm/ GW  and  0±  —  30  cm/ GW  were  obtained  for  polarisations  parallel  and  perpendic¬ 
ular  to  the  director.  In  the  isotropic  phase,  the  nonlinear  absorption  value,  0i  =  103 
cm /GW  was  obtained.  Compared  to  the  pure  material,  the  doping  increased  the  largest 
nonlinearity  by  a  factor  of  six. 

We  measured  the  OPL  characteristics  of  the  two  samples  under  simulated  test  bed 
conditions.  Both  performed  well.  Figs.  23  and  24  show  the  power  Smiting  curves  for 
these  samples. 

Increasing  the  thickness  of  the  sample  and  increasing  the  amount  of  dye  decreased  the 
transmission  to  80%  of  the  undyed  level.  A  relatively  small  improvement  is  accordingly 
observed  in  the  power  limiting  characteristics  of  sample  2.  Notice  that  the  horiaontal 
axis  is  much  enlarged  over  the  vertical  so  that  the  turn  over  of  the  transmission  curve  is 
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Figure  23:  Optical  Power  Limiting  curve  for  Sample  1 
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Figure  24:  Optical  Power  limiting  curve  for  Sample  2 

magnified.  From  this  curve  one  sees  an  onset  of  power  limiting  at  approximately  3  ftJ  and 
the  maximum  transmission  is  accordingly  approximately  3  pJ.  This  is  almost  ideal  OPL 
behavior.  It  further  illustrates  the  extent  to  which  improvements  have  progressed  towards 
operational  devices.  When  the  project  started,  very  few  materials  looked  like  they  would 
within  an  order  of  magnitude  of  this  degree  of  performance. 

Damage  in  these  samples  occurred  at  40  pJ  and  was  caused  by  the  polyimid  alignment 
layer. 


5  Picosecond  Nonlinear  Optical  Response 

The  third  time  scale  investigated  for  optica]  nonlinearities  in  liquid  crystals  was  the  pi- 
cosecond  regime.  Using  the  Nd  YAG  with  a  33  picosecond  pulse,  s-scan  measurements 
were  carried  out  to  determine  the  nonlinear  index  of  refraction  and  nonlinear  absorption  of 
liquid  crystal  materials.  Tabulated  results  covering  all  three  time  scales  are  contained  in 
Appendix  A.  Here  we  report  on  one  compound,  5CB,  which  was  studied  most  extensively. 
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We  were  able  to  obtain  eery  good  data  in  tbe  picosecond  regime.  We  attach  a  s-acan 
produced  at  the  CREOL  facility.  Thi*  is  a  email  aperture  scan  and  shows  a  mack  mailer 
venation  in  intensity  (4%  )  than  tbe  nanosecond  scans  (40%  ).  Still,  as  can  be  seen,  tbe 
scatter  in  the  data  is  very  small  and  hence  we  are  confident  in  out  measured  values  of 
nonlinear  index. 
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Figure  26:  2-scan  of  sample  with  2.3  pJ  pulses  of  33  picosecond  duration. 
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5.1  Mechanism  Analysis 


In  thin  section,  we  bring  together  the  results  of  oar  measurements  on  8  different  liquid 
crystals  and  discuss  the  conclusions  we  have  readied  about  mechanisms.  Five  of  these 
materials  are  commercially  available;  however,  OPL-7-1,  OPL-10-1  and  OPL-10-2,  were 
synthesised  by  the  group  of  L.C.  Chien  at  the  Liquid  Crystal  Institute  at  Kent  State 
University.  In  addition  to  pure  liquid  crystals,  a  dyed  sample  (anthraquinone  dye  drr 
from  BDH)  of  8CB  was  also  studied.  The  experimental  results  for  these  materials  is 
provided  in  Appendix  A,  but  we  extract  those  relevant  to  the  current  discussion  and 
present  them  here  also. 

The  molecular  structures  are  also  shown  along  with  the  wavelength,  A,  at  which  the 
measurements  were  carried  out,  6  the  laser  pulse  width,  and  T  the  sample  of  temperature. 
The  subscripts  ||  and  x  indicate  polarisation  parallel  and  perpendicular  to  the  director, 
and  i  denotes  the  isotropic  phase.  The  uncertainty  in  the  figures  is  approximately  10  %. 
Also,  in  Tables  1,  2  and  3  the  previous  results  on  5CB,  8CB,  E7  are  included  for  the 
purpose  of  comparison. 

Of  all  the  pure  materials  studied  on  the  nanosecond  time  scale,  5CB  shows  the  largest 
nonlinearities.  To  interpret  the  nanosecond  measurements,  it  is  useful  to  briefly  summarise 
the  results  of  millisecond  measurements  given  in  Table  2. 


5.1.1  Millisecond  Mechanisms 


In  our  geometry,  where  reorientation  is  not  expected  to  take  place,  the  nonlinearity  on 
the  millisecond  time  scale  originates  from  laser  heating  due  to  linear  absorption.  The 
resulting  temperature  increase  causes  a  decrease  in  the  degree  of  orientational  order  and 
in  the  density,  with  the  result  that,  in  most  materials,  the  extraordinary  index  decreases 
while  the  ordinary  index  increases.  This  gives  rise  to  n^i  <  0  and  »aX  >  0,  as  observed 
in  the  case  of  5CB.  Furthermore,— »j||,  njX  and  the  nonlinear  birefringence,  »2X  —  w^i  all 
increase  dramatically  as  the  nematic-isotropic  temperature  is  approached  from  below.  At 
the  intensities  used  in  these  measurements,  no  nonlinear  absorption  is  observed.  The  signs 
of  the  two  nonlinear  refractive  indices  measured  for  5CB  and  ZLI-2303  are  in  accordance 
with  the  above  trend.  However,  for  Z LI- 1538  and  T15  both  nonlinear  refractive  indices  are 
negative.  ZLI-1538  is  the  only  material  studied  whose  core  consists  of  cyclohexane  rather 
than  benaene  rings,  and  hence  is  without  conjugation.  Due  to  the  lack  of  conjugation, 
the  linear  polarisability  and  the  refractive  indices  are  small,  and  hence  order  parameter 
changes  are  not  effective  in  changing  the  index  of  refraction.  Thus,  the  effect  of  thermal 
expansion  or  other  smaller  nonlinearity  may  be  dominant.  T15  on  the  other  hand  has 
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Material 

m(esu) 

1 

I 

ExiSH 

5CB 

4-cyano-4'-n-pentylbiphenyl 

1:  -54x10-" 

1  265 

24 

X:  +8.3x10-" 

i:  -24x10"" 

X:  36 

i:  114 

24 

40 

8CB 

4-  cy ano-4'-  n-octy lbiphenyl 

U:  -26x10-" 

«:  246 

22 

h”c'-©-0-cn 

1:  +3.7x10-" 

X:  20 

24 

d27  doped(<0.2%)  8CB 

H:  -152x10-" 

|:  380 

22 

X:  +10.5x10-" 

X:  30 

22 

i:  -48x10-" 

i:  193 

40 

OPL-7-1 


4-(l-octynyl)-4-cyanobiphenyl 

»: 

-27x10"" 

|| :  112 

37 

X: 

• 

i : 

+  7.0x10-" 

-5.5x10-" 

X:  19 

i:  33 

37 

60 

OPL-10-1 

4-4'-dipentylazoxy  benzene 

0 

II : 

-9.9x10-" 

II :  M 

40 

X: 

i : 

+4.5x10-" 

-2.4x10-" 

X :  3.7 

i:  11 

40 

80 

OPL-10-2 

4-4 -dihexylazoxy  benzene 
n 

H: 

-15x10-" 

II :  38 

50 

C,HU 

X: 

• 

l : 

+  2.0x10"" 

+  5.0x10-" 

X:  4.2 

i:  21 

50 

66 

Table  1:  Summary  of  measurement*  of  nanosecond  nonfinearities  at  532  nm 
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Table  I  (continued) 


Material 

n*(eso) 

fi(an/GW) 

TfC) 

T15 

4-cyaao-4" -n-alkyl-p-terphenyl 

-52x10’“ 

|:  270 

180 

h"ci_ <Oh-0>—0^-cn 

X:  +21x10'“ 

JL:  18 

180 

E7 

-38x10““ 

H:  284 

24 

X:  +3.3x10““ 

X:  40 

24 

CB15 

i:  +7.4x10““ 

i:  38 

24 

ZLI-1538 

Cydohexylcydohexane 

|| :  -0 

II :  ~0  « 

70 

h^--CMI)~cn 

JL:  ~0 

i :  ~0 

X :  ~  0 

i:  -0 

70 

90 

ZLI-2303 

Phenylpyrimidine 

|| :  -20x10"“ 

II :  67 

45 

X:  +3.0x10"“ 

i:  -3.8x10”“ 

X:  8.2 

i:  19 

45 

85 

Table  1:  Continued.  Measurements  of  nanosecond  nonlinearities  at  533  nm 
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Material 

j?(cm/GW) 

TfC) 

5CB 

4-cyano-4'-n-pentylbiphenyl 

|j:  -10x10-“ 

II:  0 

24 

X:  +2.0x10”“ 

i:  -2.0x10”“ 

X:  0 

i:  a 

24 

40 

T15 

4-cyano-4" -n-alkyl-p-terphenyl 

H:  -0.68x10-“ 

»:  0 

180 

h"c,-hO— 

X:  -0.092x10-“ 

X:  0 

180 

CB15 

i:  -0.25x10"“ 

i:  0 

22 

HjC2HCH2C— ^^“CN 

ZLI-1538 

Cyclohexylcyclohexane 

H:  -0.43x10-“ 

II  ••  « 

70 

H*-Q - 

X:  -0.13x10"“ 

i:  -0.35x10”“ 

X:  0 

i:  0 

70 

97 

ZLI-2303 

Phenylpyrimidine 

|| :  -0.97x10-“ 

H  ••  o 

45 

huc« — — ^3) — °c<Hi3 

X:  +0.082x10-“ 

i:  -0.3x10*“ 

X:  0 

i:  0 

45 

85 

Table  2:  Summary  of  measurements  of  nonlinearities  at  514  nm,  10  ms  pulse 
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Material 

na(«u) 

fi(aa/GVf) 

TfC) 

5CB 

4-cyano-V-n-pentylbiphenyl 

II :  +  1.04xl0"u 

|| :  2.26 

24 

h,^hOmO>- 

J.:  +0.69x10” 11 

-L:  0.78 

24 

Table  3:  Summary  of  measurements  of  picosecond  nonlinearities  at  532  nm 

three  bensene  rings  instead  of  two,  and  thus  the  reason  why  both  nonlinear  indices  are 
negative  is  not  clear.  The  observed  nonlinearities  seem  to  be  relatively  insensitive  to 
structural  changes  other  than  the  presence  of  bensene  rings.  For  example,  the  addition  of 
a  chiral  branched  alkyl  chain  in  CB15  does  not  change  the  nonlinear  index  in  the  isotropic 
phase  appreciably. 


5.1.2  Nanosecond  Mechanisms 

Next,  we  consider  possible  mechanisms  responsible  for  the  nonlinearities  on  the  nano¬ 
second  time  scale.  Strong  nonlinear  refraction,  nonlinear  birefringence  and  nonlinear 
absorption  is  observed  in  these  materials.  By  noting  that  the  temperature  dependence  of 
the  nonlinear  indices  on  the  nanosecond  time  scale  is  fundamentally  different  from  that 
on  the  millisecond  scale,  we  conclude  that  the  dominant  mechanism  is  not  laser  heating 
due  to  linear  absorption.  This  is  further  substantiated  by  the  observations  that  m  for 
GB15  and  nx  for  T15  change  sign  as  the  pulse  width  is  reduced  from  ms  to  ns. 

Although  many  other  mechanisms,  such  as  direct  optical  field  induced  orientational  ef¬ 
fects,  electrostrictive  effects,  elect  restrictive  effects  combined  with  shear  alignment,  may 
contribute  to  the  observed  nonlinearities,  order  of  magnitude  calculations  based  on  simple 
theoretical  models  indicate  that  they  can  not  be  the  primary  cause.  The  intensity  depen¬ 
dent  measurements  on  5CB  indicate  that  n*  is  independent  of  intensity,  and  therefore  it 
is  a  third  order  process.  Analysis  of  the  measured  peak  and  valley  positions  on  the  Z-scaa 
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curvet  are  consistent  with  this.  Two-pulse  measurements  indicate  that  the  mechanism 
is  slow,  that  is,  the  response  time  can  not  be  less  than  about  5ns.  One  possibility  is 
that  the  nonlinear  refraction  originates  from  an  excited  state  associated  with  the  linear 
one-photon  absorption.  On  the  other  hand,  for  5CB,  intensity  dependence  measurements 
show  that  the  nonlinear  absorption  coefficient,  is  nearly  proportional  to  the  intensity. 
This  is  a  signature  of  the  5th-order  process;  two  pulse  measurements  furthermore  show 
that  the  associated  response  time  is  less  than  about  5ns.  It  appears  likely  therefore,  that 
at  least  in  the  case  of  5CB,  two  different  processes  are  operating  simultaneously,  a  slow 
third-order  process  which  contributes  primarily  to  the  nonlinear  refraction,  and  a  faster 
5th-order  process  which  contributes  primarily  to  the  nonlinear  absorption. 

The  nonlinear  behavior  is  linked  to  the  presence  of  the  benzene  rings  in  the  molecules. 
For  example,  for  ZLI-1538  no  observable  nonlinearities  (either  nonlinear  refraction  or 
absorption)  have  been  detected;  this  material  does  not  have  a  benzene  ring.  Conversely, 
the  addition  of  the  anthraquinone  dye  to  8CB  significantly  enhances  the  observed 
nonlinearity,  without  changing  its  character  (i.e.  signs  of  nj  are  unchanged  from  that  of 
the  pure  material;  both  nonlinear  refraction  and  nonlinear  absorption  increase). 

Several  other  features  of  molecular  structure  of  our  selection  of  liquid  crystals  did  not 
play  an  important  role  in  determining  the  size  of  the  nonlinearity.  As  can  be  seen  in  Table 
1,  the  addition  of  the  acetylene  group  (OPL-7-1)  does  not  significantly  contribute  to  the 
nonlinearity.  Similarly,  varying  the  linkage  by  adding  small  conjugated  structures  (azoxy 
group)  also  has  little  effect  on  the  nanosecond  response. 

Our  key  findings  on  nanosecond  time  scale  experiment  are  the  following:  the  nonlinear 
refraction  in  a  number  of  liquid  crystals  in  anomalously  large,  the  nonlinear  refraction 
is  accompanied  by  large  nonlinear  absorption,  and  the  magnitude  of  the  nonlinearity  is 
closely  linked  to  the  number  of  benzene  rings  in  the  molecules,  and  appears  insensitive 
to  the  presence  of  other  conjugated  units.  Our  results  indicate  that  in  5CB  a  slow  third- 
order  process  is  responsible  for  the  nonlinear  refraction,  while  a  faster  fifth-order  process 
is  responsible  for  the  nonlinear  absorption. 


5.1.3  Picosecond  Mechanisms 


If  the  measured  0  was  due  to  three  photon  absorption,  then  picosecond  measurements 
should  show  very  large  values  of  /?,  since  the  intensities  are  more  than  one  order  of 
magnitude  higher.  However,  0  measured  with  33ps  is  smaller;  values  are  given  in  Table 
3.  It  is  likely  therefore  that  the  response  time  of  the  absorption  process  is  long  compared 
to  picoseconds  (but  short  compared  to  10ns)  and  this  accounts  for  the  decrease  of  0 
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ceen  in  the  picosecond  measurements.  Such  a  process  could  be  excited  state  absorption 
from  a  two-photon  excited  state.  Since  the  nonlinear  refraction  which  dominates  on  the 
nanosecond  time  scale  is  slow,  it  will  be  attenuated  by  more  than  two  orders  of  magnitude 
in  the  picosecond  measurements.  The  picosecond  n3  values  thus  originate  from  another 
mechanism,  such  as  nonresonant  electronic  hyperpolarizability.  The  result  that  both  n3 
and  are  very  small  at  A  =  1.06pm  is  not  consistent  with  this  model.  The  contribution  of 
the  excited  state  at  this  wavelength  to  the  molecular  polarizability  could  be  much  smaller 
than  that  at  532  nm,  and  the  multiphoton  absorption  cross-section  is  certain  to  be  very 
different. 


5.1.4  Overall  Mechanism  Discussion 

There  are  a  variety  of  other  possible  origins  to  nonlinear  response  in  liquid  crystals  in  the 
nanosecond  regime.  We  have  considered  above  those  most  likely  to  be  large  and  those 
which  are  clearly  evident  from  the  experimental  measurements.  Still,  we  are  convinced 
that  other  nonlinearities  are  present  and  their  understanding  and  quantification  are  nec¬ 
essary  to  actually  predict  theoretically  the  nonlinear  behavior  over  a  wide  range  of  time 
scales,  wavelengths  and  device  geometries. 

In  particular,  in  the  nanosecond  regime,  it  is  possible  that  electrostriction,  short-axis 
reorientation,  photostimulated  conformational  changes  and  complex  intramolecular  en¬ 
ergy  transfers  play  a  role,  but  more  study  is  needed  to  understand  this  behavior.  It  seems 
to  be  true  that  while  the  nanosecond  OPL  properties  of  several  liquid  crystals  are  large 
and  operative  over  many  wavelengths  on  the  nanosecond  regime,  this  cannot  be  attributed 
to  a  «ingl*  third  order  nonlinear  process  and  we  do  not  understand  the  factors  that  make 
the  materials  work  well. 

Even  in  the  millisecond  regime  where  the  mechanisms  involving  thermally  induced  dis¬ 
order  and  optical  field  reorientation  seem  to  be  relatively  clear,  specifics  for  each  molecule 
are  however  not  quantified. 

Only  the  picosecond  regime  seems  to  be  simple  to  understand,  but  it  remains  difficult 
to  predict  for  any  given  molecule.  We  feel  secure  that  the  picosecond  nonlinearities  are 
almost  entirely  electronic  but,  while  we  have  a  rule  of  thumb  ability  to  predict,  there 
remains  an  inability  to  predict  accurately  and  there  remains  a  number  of  molecules  which 
have  not  been  measured. 
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6  Damage  Studies 


In  both  the  CW  and  nanosecond  OPL  studies,  damage  to  the  samples  was  observed  at  high 
intensity.  Studies  were  undertaken  to  identify  the  damage  thresholds  and  mechanisms. 


6.1  CW  Damage  Studies 

The  damage  threshold  studies  were  motivated  by  the  idea  that  there  may  be  some  focusing 
of  the  laser  beam  by  the  liquid  crystal  droplets  causing  a  damage  threshold  in  the  PDLC 
to  be  lower  than  that  inherent  for  either  the  polymer  host  or  liquid  crystal  materials 
separately.  Controlled  experiments  were  thus  performed  on  the  three  types  of  sample, 
each  10  pm  thick,  sandwiched  between  ITO  coated  glass  plates. 

The  sample  types  were: 


•  Type  I  PDLC  (Polymer  with  33%  liquid  crystal  by  weight) 

•  Type  II  Pure  Polymer  (Epon  plus  Capcure  in  equal  parts) 

•  Type  III  Pure  liquid  crystal  (E7,  commercial  mixture) 


The  following  thresholds  were  observed: 


•  Type  I:  26±  lkW/cm2  (10ms  pulse,  95  pm  beam  diameter) 

•  Type  I:  18±  lkW/cm2  (20ms  pulse,  95  pm  beam  diameter) 

•  Type  II:  32±  lkW/cm2  (10ms  pulse,  95pm  beam  diameter) 

•  Type  III:  >  40kW/cm2  (10ms  pulse,  95  pm  beam  diameter) 

It  is  clear  that  the  liquid  crystal  material  is  not  the  damaged  material,  but  the  fact 
that  PDLC  damages  at  a  lower  threshold  than  the  polymer  suggests  that  some  mechanism 
such  as  focusing  is  operational.  The  reduction  in  threshold  with  the  addition  of  the  liquid 
crystal  is  a  small  effect,  20%,  but  the  mechanism  is  not  completely  understood.  The 
criterion  applied  to  determine  damage  was  discoloration  (brown).  Both  PDLC  and  pure 
polymer  samples  showed  the  same  type  of  damage. 


Following  an  informal  report  on  results  of  measurements  by  David  Pepper  at  Hughes, 
we  investigated  the  role  played  by  the  ITO  coating.  Similar  samples  (thickness  20pm)  of 
the  three  types  were  studied  between  uncoated  glass.  The  results  were: 


•  Type  I:  >  41kW/cm2  (10ms  pulse,  95pm  beam  diameter) 

•  Type  II:  >  41kW / cm2  (10ms  pulse,  95pm  beam  diameter) 

•  Type  III:  >  41kW/cm2  (10ms  pulse,  95  pm  beam  diameter) 


The  6%  optical  absorption  by  the  ITO  is  apparently  large  enough  to  contribute  the  pri¬ 
mary  effect  in  damaging  the  polymer  material.  This  absorption  is  essentially  unnoticeable 
in  normal  electro-optic  applications  involving  liquid  crystals. 

These  experiments  show  that  the  ITO  is  the  vehicle  whereby  damage  is  inflicted,  but 
the  difference  between  PDLC  threshold  and  polymer  threshold  remains  unexplained,  but 
is  likely  due  to  the  difference  in  heat  diffusion  away  from  the  surface. 


6.2  Nanosecond  Laser  Damage 

We  have  studied  laser  damage  in  homeotropidy  (parallel)  aligned  samples  of  5CB.  We 
found  that  the  cells  are  damaged  by  nanosecond  laser  pulses  at  relatively  low  fluence  (~1 
J/cm2)  due  to  self-defocusing  of  the  beam  by  the  fiquid  crystal  sample.  The  damage  is 
produced  on  the  interior  surface  of  the  glass;  there  is  no  apparent  damage  to  the  liquid 
crystal.  This  effect  may  be  utilised  as  a  sacrificial  optical  fuse. 

With  the  materials  deliver  for  the  ARDEO  DVO  test  bed,  we  found  that  damage  in 
these  samples  occurred  at  40  pj  and  was  caused  by  the  polyimid  alignment  layer. 

We  can  not  confirm  our  final  evaluation  of  the  damage  mechanisms,  but  we  have  several 
thoughts.  It  appears  that  in  the  nanosecond  regime,  the  damage  condition  is  simply  that 
they  intensity  at  the  interface  of  the  sample  cell  must  reach  the  damage  threshold  of  the 
material  of  the  surface.  In  the  case  when  there  is  an  alignment  layer  such  as  polyimid, 
the  polyimid  is  damaged.  If  the  interface  is  simply  a  liquid  crystal/glass  interface,  then 
the  glass  is  damaged.  We  have  no  independent  determination  of  the  damage  threshold  of 
polyimid,  but  for  glass,  it  appears  that  the  damage  intensity  is  simply  the  usual  damage 
level  of  a  glass  surface  and  has  nothing  to  do  with  the  liquid  crystal.  The  only  effect 
of  the  liquid  crystal  seems  to  be  that  the  focal  point  or  location  of  the  beam  waist  is 
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modified  as  a  function  of  intensity  due  to  the  presence  of  the  liquid  crystal.  If  the  beam 
waist  is  defocused  from  the  center  of  the  sample  to  the  surface,  then  the  surface  damage 
resuts  sooner  than  it  does  in  a  thicker  sample  or  a  sample  without  liquid  crystal.  This 
can  obviously  be  incorporated  in  the  design  of  a  device  by  adjusting  the  device  geometry 
(thickness)  and  choosing  high  threshold  materials  for  the  sample  holder. 
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7  Theory  and  Modeling 


The  modeling  and  theory  effort  under  this  project  undertook  to  provide  predictive  support 
for  materials  selection  and  theoretical  modeling  to  explain  the  mechanisms  involved  in 
the  nonlinear  response  of  liquid  crystals  materials.  Three  separate  areas  of  activities 
were  involved,  electronic  structure  modeling,  optical  fields  induced  order  modeling,  and 
nonlinear  optical  propagation  modeling. 


7.1  Electronic  Structure  Modeling 

While  there  is  no  set  of  rules  which  specify  which  molecules  will  actually  exhibit  liquid 
crystalline  phases,  there  are  strong  chemical  trends  which  are  well  known  and  extensively 
studied.  Molecular  polarizability  is  an  essential  quality  of  most  liquid  crystals.  All  but  a 
few  liquid  crystalline  molecules  have  an  aromatic  core  usually  composed  of  one  or  more 
benzene  rings.  In  addition,  it  is  well  known  that  the  liquid  crystalline  properties  are 
strongly  influenced  by  substituent  groups  attached  to  this  core. 

Molecular  polarizability  strongly  correlates  with  electronic  nonlinear  optical  properties 
of  organics.  That  is  why  one  should  expect  liquid  crystals  to  exhibit  strong  nonlinear 
optical  properties.  Given  the  variety  of  liquid  crystals,  one  should  also  expect  a  variety, 
and  some  of  them  large,  nonlinear  optical  characteristics.  These  qualitative  considerations 
are  reasons  to  study  nonlinear  optical  properties  of  liquid  crystals,  but  they  do  not  provide 
an  explanation  of  those  properties.  Ultimately,  as  we  have  discovered,  there  are  properties 
of  the  liquid  crystal  phases  which  affect  the  nonlinear  properties  and  there  are  a  number 
of  mechanisms  for  optical  nonlinearities  that  interact  in  a  complex  manner. 

While  molecular  polarizability  correlates  with  nonlinear  properties,  it  is  not  the  origin 
of  the  nonlinearity.  It  is  the  electronic  hypcrpolarizability  that  is  the  origin  of  molecular 
electronic  nonlinearities.  This  requires  a  quantum  mechanical  theory  of  molecular  elec¬ 
tronic  structure.  As  our  first  theoretical  effort  was  to  investigate  electronic  nonlinearities, 
we  adapted  a  fairly  direct  method  of  electronic  structure  calculation  to  use  to  determine 
trends  in  and  origins  of  electronic  nonlinearities  in  liquid  crystals.  The  theoretical  ap¬ 
proach  was  to  isolate  the  electronic  structure  of  the  aromatic  core  of  the  liquid  crystal 
molecule,  since  past  experience  has  shown  us  that  this  molecular  segment  contribute  dom¬ 
inantly  to  the  linear  and  nonlinear  molecular  polarizability.  We  employed  the  Huckel  and 
PPP  models  to  describe  the  electronic  structure  of  the  core.  This  models  includes  only 
the  so-called  pi  electrons,  i.e.  those  involves  in  the  conjugation  of  double  and  triple  bonds 
in  the  molecule.  Substituent  groups  and  heteroatoms  are  treated  as  perturbation  of  the 
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banc  Hamiltonian. 

Additional  background  on  the  specific*  of  the  quantum  mechanical  calculation  for  this 
portion  of  the  theoretical  effort  can  be  found  in  the  Quarterly  Scientific  and  Technical 
Reports.  Also  in  the  Report  are  extensive  listings  of  the  many  molecules  calculated  and 
discussions  of  trends  found.  Here,  we  will  summarise  a  few  aspects  that  were  used  in 
guiding  materials  selection  and  synthesis. 

7.1.1  Small  Molecule  Comparisons 

We  will  first  discuss  our  efforts  at  verification  of  the  electronic  structure  work.  We  did 
this  by  comparing  the  results  of  calculations  of  the  two  different  Hamiltonians.  The 
Huckel  Hamiltonian  is  quite  simple  and  amenable  to  calculation,  even  for  large  molecules. 
Its  simplicity  makes  it  suspect  however.  The  PPP  Hamiltonian  is  more  complex  and 
includes  electron-electron  interactions.  It  is  computationally  difficult  because  the  cost  of 
a  calculation  increases  as  the  factorial  of  the  number  of  electrons.  We  reasoned  that  if 
the  complicated  Hamiltonian  predict  the  same  trends  in  nonlinear  properties  as  does  the 
simple  Hamiltonian,  then  we  could  rely  on  the  predictions  of  the  Huckel  Hamiltonian  in 
developing  materials  selection  criteria. 

The  largest  molecules  where  we  could  use  both  Hamiltonians  were  substituted  benzene 
molecules  and  straight  chain  alkenes.  Fig.  26  shows  the  results  of  a  series  of  calculations 
of  the  molecular  hyperpolarizability  of  substituted  benzenes.  Included  are  both  sets  of 
theoretical  results  and  the  experimental  results  taken  from  the  literature. 

The  agreement  in  Fig.  26  is  quite  gratifying.  Not  only  do  the  two  theoretical  predic¬ 
tions  have  the  same  trend  as  each  other  and  the  experimental  results,  the  quantitative 
agreement  is  reasonably  good  among  all  three.  The  extend  to  which  there  is  disagree¬ 
ment  can  be  understood  semiquantitatively  in  terms  if  the  hyperpolarizability  of  the  S 
electrons  and  sigma  bonds  that  are  not  included  in  the  theories.  Such  a  discussion  has 
been  published [20],  but  since  we  are  interested  in  trends  and  large  changes  in  hyperpo¬ 
larizability,  we  will  not  include  that  discussion  here.  These  considerations  were  discussed 
extensively  in  the  Quarterly  Scientific  and  Technical  Reports. 
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Halogen  substitution 

Figure  26:  Hyperpolarizability  of  substituted  benzenes 
7.1.2  Lurge  Molecule  Cnlculsitions 

We  calculated  the  linear  (first  order)  polarizabilities  and  the  second  and  third  order  elec¬ 
tronic  hyperpolarizabilities  for  a  number  of  liquid  crystal  molecules.  Calculations  were 
performed  for  a  good  many  small  molecules  as  well,  as  part  of  the  verification  process. 

Fig.  27  shows  the  chemical  structure  of  the  aromatic  cores  for  a  number  of  liquid 
crystal  molecules  that  we  considered.  We  provide  this  as  a  guide  for  the  tables  which 
follow. 
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Figure  27:  Structure  end  names  of  molecules,  used  in  tables 
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biphenyl 

9. SOM 

2.9707 

1.9601 

stllbene 

5.4219 

2.9519 

2.8581 

azobenzene 

5.3200 

2.9747 

2.1322 

diphenylacecylene 

7.2649 

2.9804 

3.2220 

cyanoblphenyl 

8.6129 

2.9609 

3.7245 

terphenyl 

6.9794 

9.5360 

3.9097 

azoxybenzene 

6.5519 

9.9954 

3.5160 

diphenyl-butadiene 

8.5798 

4.3919 

4.2672 

dipheny 1 -butediyne 

12.963 

2.9812 

5.1221 

Table  4:  Molecular  linear  polarisabilities 

Af  indicated  in  the  introductory  discussion,  linear  polarizability  is  expected  to  follow 
nonlinear  properties  for  several  reasons.  The  nonlinear  properties  in  the  picosecond  repine 
are  expected  to  be  electronic  and  largely  caused  by  the  hyperpolarisabiHty  of  the  aromatic 
core.  As  we  will  see,  the  size  of  the  core  is  a  major  factor  in  both  the  linear  and  hyper 
polarizabilities,  hence  the  correlations.  We  refer  here  to  nonruonant  nonlinearities.  The 
presence  of  fortuitously  located  excited  states  can  cause  enormous  nonlinearities  over  a 
very  narrow  wavelength  range.  Our  charge  was  to  seek  materials  suitable  for  "frequency 
agile”  lasers,  thus  resonant  phenomena  were  not  emphasised.  The  polarisabilities  reported 
here  are  static,  hence  we  believe  more  appropriate  for  predicting  trends  in  materials  over 
a  wide  wavelength  regime. 

Table  4  shows  the  linear  polarizabilities  calculated  for  the  molecules  under  discussion. 
It  is  followed  by  Table  5  showing  the  third  order  hyperpolarizabilities  calculated. 

The  table  on  molecular  linear  polarisabilities  gives  the  values  of  the  nonzero  compo¬ 
nents  of  the  molecular  polarizability  tensor  a.  These  values  are  due  only  to  the  conjugated 
(aromatic)  portion  of  the  molecule.  The  long-axis  polarizabilities  track  with  the  nematic 
polarization  in  the  direction  of  the  molecules  and  the  short  axis  polarizability  tracks  with 
the  perpendicular  direction  in  the  molecule.  We  will  ignore  of  long  ™  polarisa¬ 

tion  with  short  axis  polarisation  in  the  nematic  phase  caused  by  thermal  averaging,  but 
in  the  isotropic  phase  this  averaging  is  complete  and  it  is  the  average  that  tracks  with  the 
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Table  5:  Molecular  hyperpolarisabilities 


isotropic  phase  polarisability  (i.e.  index  of  refraction). 

The  Table  on  molecular  hyperpolarisabilities  gives  the  values  of  the  nonaezo  compo¬ 
nents  of  the  molecular  hyperpolarisabaKty  tensor  7.  The  hyperpolarisability  is  a  more 
complex  entity  than  the  linear  polarisability  once  it  is  a  rank  four  tensor  (rather  than 
rank  two).  These  values  are  also  due  only  to  the  conjugated  (aromatic)  portion  of  the 
molecule.  The  long-axis  hyperpolarisabilities  track  with  the  nematic  hyperpolarisation 
in  the  direction  of  the  molecules  and  the  short  axis  hyperpolarisability  tracks  with  the 
perpendicular  direction  in  the  molecule.  We  will  ignore  ™«i«g  of  long  axis  hyperpolar¬ 
isation  with  short  axis  polarisation  in  the  nematic  phase  caused  by  thermal  averaging, 
but  in  the  isotropic  phase  this  averaging  is  again  needed  but  is  more  complex  that  the 
linear  case.  The  entry  for  the  average  has  been  calculated  and  it  tracks  with  the  isotropic 
phase  hyperpolarisability  (i.e.  nonlinear  index  of  refraction). 


7.1.3  Discussion  of  Electronic  Nonlinearities 


One  can  see  the  origin  of  our  selection  rules  in  the  tabulated  results.  The  experiments 
in  the  millisecond  regime  (CW)  will  be  most  influenced  by  the  linear  polarisability.  This 
is  because  the  induced  dipole  is  directly  proportional  to  the  optical  torque,  and  that 
portion  of  the  dipole  which  is  effective  in  optical  reorientation  is  directly  proportional  to 
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ike  anisotropy  of  the  linear  polariaahflity  tensor,  Aa.  The  net  effective  nonfinear  index 
resulting  from  optical  reorientation  should  then  track  with  the  anisotropy. 

Simply  by  looking  at  the  molecules,  it  can  be  seen  that  the  anisotropy  in  polarisability 
tracks  with  the  anisotropy  in  the  shape  of  the  molecular  core.  This  rule  of  thumb  is  easily 
implemented  as  a  synthesis  guide.  This  particular  rule  has  advantages  and  disadvantages. 
On  the  positive  side,  the  tendency  of  a  molecule  to  form  a  Hquid  crystal  is  well  known 
to  increase  as  the  anisotropy  in  the  core  increases.  On  the  negative  side,  there  is  also 
generally  an  increase  in  the  temperature  at  which  the  nematic  forms,  sometime  getting 
well  above  room  temperature,  and  there  is  an  increase  in  viscosity  making  it  less  liquid 
and  more  difficult  to  reorient. 

Fortunately,  much  of  the  above  discussion  applies  to  the  general  trends  observed  for 
hyperpolarisability.  Indeed,  hyperpolarisability  increased  even  more  rapidly  than  linear 
polarisability  when  compared  to  the  the  geometric  anisotropy  of  the  aromatic  core.  The 
more  complex  nature  of  the  hyperpolarisability  make  such  rules  of  thumb  less  clearly 
applicable.  In  particular,  the  off  diagonal  elements  of  the  hyperpolarisability  are  often 
large  and  of  opposite  sign  to  the  long  axis  hyperpolarisability.  The  assumption  of  ignoring 
thermal  averaging  in  the  nematic  phase  is  thus  compromising.  While  it  is  not  difficult 
to  carry  out  thermal  averaging  (and  we  have  done  this),  it  is  not  reliable  because  the 
absolute  magnitudes  of  the  components  of  the  hyperpolarisability  tensors  are  not  reliable. 
As  we  stated  in  the  beginning,  we  are  seeking  trend#  from  the  Huckel  model,  not  absolute 
prediction,  when  we  use  it  for  large  molecules.  In  particular,  we  believe  that  the  Huckel 
model  significantly  over  estimates  7XSW  and  it  enters  in  a  crucial  way  in  the  averaging 
scheme. 

Rather  than  use  the  averages  that  we  can  extract  from  the  hyperpolarisability  tensors, 
we  adopted  a  more  qualitative  approach.  One  can  see,  for  example,  that  the  addition  of 
double  bonds  between  bensene  rings  increases  some  hyperpolarisability  components.  It 
does  not  however  increase  the  dominant  component  as  much  as  the  addition  of  another 
bensene  ring.  This,  simply  looking  at  the  shape  of  the  molecule  is  not  nearly  as  adequate 
in  predicting  hyperpolarisability  as  loving  at  the  actual  calculations.  These  ideas  were 
to  be  testing  experimentally,  and  that  was  part  of  the  reason  for  choosing  to  synthesise 
some  materials  with  double  bonds  connecting  or  appended  to  bensene  ring  systems. 

These  calculations  told  us  that  by  varying  structure  of  the  core  and  varying  linkages 
between  cores  and  by  varying  substituent  groups,  we  could  expect  to  see  variations  in 
the  nonlinear  properties  of  the  liquid  crystals.  While  our  discussion  here  is  not  complete, 
and  it  was  supplemented  with  heuristic  expectations  of  molecular  origins  of  nonlinear 
response,  it  did  help  us  to  choose  to  select  materials  for  evaluation.  We  extract  here  a 
portion  of  the  discussion  on  Appendix  A  where  we  state  the  essential  elements  of  our 
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adaption  rules. 


•  Varying  uses  of  conjugated  core* 

•  Varying  linkages  between  benaene  rings 

•  Varying  substituent  groups 

•  Ease  of  synthesis  or  even  possible  commercial  availability. 

•  Simple  and  small  conjugation  structures  to  facilitate  theoretical  almUtwin 

In  Sec.  8.2  we  list  the  molecules  chosen  for  study.  Many  of  these  were  synthesised  in 
house.  To  obtain  a  representative  sample  covering  a  reasonable  portion  of  the  above  list 
we  had  to  synthesise  some  even  though  they  represented  simpler  structures  than  some 
commercially  available  compounds.  In  the  end,  we  had  mote  compounds  than  we  had 
time  to  measure  and  we  had  identified  dosens  more  (Appendix  A)  that  were  not  acquired. 
The  extension  of  a  survey  involving  these  compounds  would  be  desirable. 

Of  the  compounds  selected  for  study,  one  can  see  the  application  of  the  above  rules 
and  the  theoretical  calculations.  For  instance,  the  addition  of  a  double  bond  between  two 
bensene  rings  had  a  much  smaller  predicted  effect  than  the  addition  of  a  cyano  group  on 
the  end.  Indeed,  before  we  ever  measure  the  nonlinear  response  of  5CB,  our  calculations 
(See  Table  5)  predicted  the  the  core  of  5GB  would  have  the  largest  hyperpolarisabifity. 
Only  with  the  addition  of  three  intermediate  carbons  (and  two  double  bonds)  does  the 
molecular  anisotropy  increase  enough  to  yield  a  predicted  nematic  hyperpoiarisability  as 
large  as  that  of  5GB. 

We  do  not  want  to  make  too  much  out  of  this  agreement  between  the  measure  nonlinear 
response  and  the  theoretical  predictions  of  hyperpolarismbilty.  There  are,  as  we  have 
repeated  stated,  other  nonlinear  mechanisms  involved.  Still,  we  found  these  calculations 
useful  in  determining  the  compounds  we  would  synthesise  and  purchase. 

Finally,  we  emphasise  that  we  have  discussed  here  molecular  properties.  Modeling 
these  systems  on  the  macroscopic  scale  was  a  valuable  additional  activity  to  determine 
how  the  various  nonlinear  mechanisms  would  manifest  themselves  experimentally.  Secs. 
7.2  and  7.3  describe  our  theoretical  efforts  in  those  areas. 
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7.2  Modeling  Molecular  Orientational  Relaxation  and  the  Ther¬ 
mal  Nonlinear  Refraction  of  Nematic  Liquid  Crystals 


Laser  heating  is  known  to  make  an  important  contribution  to  the  nonlinear  refraction  of 
liquid  crystals  due  to  the  large  rate  of  change  of  the  refractive  index  with  temperature, 
dn/dT.  There  is  a  strong  enhancement  of  dn/dT  near  the  nematic-isotropic  phase  transi¬ 
tion  which  is  associated  with  the  rapid  variation  of  the  nematic  order  parameter,  5,  with 
temperature.  This  has  been  seen  to  lead  to  a  strong  enhancement  of  the  nonlinear  refrac¬ 
tion.  However,  if  the  nonlinear  refraction  is  probed  with  laser  pulses  on  the  nanosecond 
time  scale,  incomplete  relaxation  of  the  order  parameter  during  the  pulse  results  in  a 
reduction  of  the  observed  nonlinear  refraction.  The  time  scale  for  relaxation  of  the  order 
parameter  also  shows  a  large  peak  near  the  phase  transition  temperature,  so  that  for  short 
laser  pulses,  the  temperature  dependence  of  the  effective  nonlinear  refractive  index,  n2, 
will  be  strongly  influenced  by  the  relaxation  dynamics  of  the  order  parameter.  We  have 
developed  a  simple  model  for  the  time  evolution  of  the  orientational  order  of  a  nematic 
liquid  crystal  subjected  to  heating  by  a  gaussian  laser  pulse.  By  treating  the  orientational 
degrees  of  freedom  with  a  phenomenological  Landau-de  Gennes  theory,  and  assuming  that 
the  other  degrees  of  freedom  equilibrate  quickly  compared  to  the  pulse  duration,  it  is  pos¬ 
sible  to  solve  for  the  time  dependence  of  the  order  parameter,  and  thus  the  indices  of 
refraction.  By  defining  an  intensity  weighted  average  of  the  change  in  index  of  refraction, 
an  effective  nonlinear  refractive  index  can  be  extracted,  such  as  would  be  measured  in 
an  experiment  of  the  z-scan  type.  Using  parameters  which  have  been  measured  for  the 
material  5CB,  it  is  found  that  the  slow  relaxation  completely  cancels  the  enhancement  of 
dn/dT  near  the  phase  transition,  resulting  in  an  effective  n2  which  is  nearly  independent 
of  temperature,  for  pulse  durations  less  than  a  few  tens  of  nanoseconds. 


7.2.1  Time  Scale  of  Thermal  Nonlinearities 


The  nonlinear  optical  response  of  liquid  crystals  has  been  the  subject  of  much  investigation 
recently. [2] ~  [9]  In  particular,  as  part  of  our  effort,  we  have  undertaken  to  measure  the 
nonlinear  optical  properties  of  various  liquid  crystalline  materials,  under  a  variety  of 
temperatures,  optical  pulse  durations,  and  optical  field  intensities  as  described  in  Secs.  3, 
4  and  5  and  published  works[4,  5].  These  experiments  made  use  of  the  z-scan  technique 
(10, 11]  which  relies  on  the  self-refraction  of  a  single  pulse.  The  results  of  the  measurements 
allow  the  samples  to  be  characterized  by  a  nonlinear  refractive  index,  n2,  where  the 
intensity  dependent  index  of  refraction  is 

n(£?2)  =  n  +  (1) 

U 
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and  a  nonlinear  absorption  coefficient,  0,  where  the  intensity  dependent  absorption  is 
given  by 

a(I)  =  a  +  0I  (2) 

One  of  the  simplest  mechanisms  which  can  contribute  to  the  nonlinear  refraction  is 
a  thermal  mechanism,  wherein  linear  absorption  of  energy  from  the  optical  field  leads 
to  heating  of  the  sample,  thereby  changing  iti  index  of  refraction.  This  mechanism  is 
expected  to  lead  to  large  effects  in  liquid  crystals  due  to  the  rapid  variation  of  the  index 
of  refraction  with  the  temperature  in  these  materials.  This  is  especially  true  near  the 
nematic-isotropic  phase  transition,  where  the  variation  of  n  with  T  is  most  rapid.  This 
effect  has  been  found  to  contribute  significantly  to  a  number  of  measurements;  more 
detailed  references  and  brief  reviews  are  published  [2,  3). 

An  estimate  of  the  effective  value  of  ns  due  to  the  thermal  effect,  as  measured  with  a 
pulsed  laser,  can  easily  be  obtained.  For  a  sample  with  a  linear  absorption  coefficient  a, 
and  heat  capacity  per  unit  volume  C,  in  a  laser  pulse  of  peak  intensity  J0  and  duration  r, 
the  change  in  temperature  is  AT  »  aI0r/C.  The  effective  value  of  na  is  approximately 

r*Zf/  *  An/E* 

»  ~ATcn/I0 

dn  otren  ( 

”  *f~C~  ^ 

where  c  is  the  speed  of  light  and  n  is  the  linear  index  of  refraction.  All  factors  of  order 
unity  have  been  dropped  for  the  purposes  of  this  order  of  magnitude  estimate.  Eq.  3 
ignores  the  diffusion  of  heat.  In  the  experiments  under  consideration,  the  laser  pulse 
duration  is  much  shorter  than  the  time  scale  for  diffusion  of  heat  across  the  beam,  so  that 
diffusion  of  heat  will  be  neglected  in  all  the  follows.  Note  that  the  result  of  Eq.  3  depends 
explicitly  on  the  pulse  duration,  r,  as  well  as  on  material  properties.  This  is  the  reason 
for  referring  to  the  “effective*'  »j. 

It  was  concluded  in  Sec.  3.3.1  and  also  published  by  Palfiy  [2,  4],  that  the  observed 
nonlinear  refraction  for  pulses  of  several  millisecond  duration  was  due  to  the  thermal 
mechanism.  This  conclusion  was  based  on  comparison  of  the  measurements  with  order 
of  magnitude  estimates  from  Eq.  3  and,  most  important,  the  fact  that  the  measured  ns 
has  a  pronounced  dependence  on  the  temperature  as  the  nematic-isotropic  transition  is 
approached  from  below.  This  strong  temperature  dependence  is  expected,  since  dn/dT  is 
strongly  temperature  dependent  for  liquid  crystals  near  the  transition  temperature.  On 
the  other  hand,  it  was  found[5]  that  for  pulses  of  7  ns  duration,  the  observed  ns  was 
essentially  independent  of  temperature,  right  up  to  the  transition.  This  was  taken  as  an 
indication  that  the  dominant  effect  on  that  time  scale  did  not  have  a  thermal  origin. 
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It  should  be  pointed  out,  however,  that  Eq.  3  is  only  appropriate  if  the  heating  is  quasi¬ 
static,  that  is  if  the  temperature  changes  slowly  enough  that  the  instantaneous  change  in 
index  of  refraction  is  related  to  the  instantaneous  change  in  temperature  by  the  equilib¬ 
rium  coefficient,  jjL  This  will  only  be  the  case  if  the  relaxation  process  which  governs 
the  change  in  index  of  refraction  takes  place  on  a  time  scale  which  is  short  compared  to 
the  pulse  duration.  If  slow  relaxation  prevents  the  index  of  refraction  from  keeping  pace 
with  the  instantaneous  temperature  during  the  optical  pulse,  then  the  observed  n2  will  be 
reduced  from  the  quasi-static  prediction.  For  nematic  liquid  crystals,  the  variation  of  the 
index  of  refraction  with  temperature  is  primarily  due  to  the  change  of  the  orientational 
order  parameter,  5,  which  is  a  measure  of  the  degree  of  alignment  of  the  long  molecular 
axes,  and  is  proportional  to  the  dielectric  anisotropy.  It  is  the  rapid  variation  of  the  order 
parameter  with  temperature  near  the  transition  which  gives  rise  to  the  large  value  of 
dn/dT.  In  a  sample  subjected  to  rapid  heating  by  a  laser  pulse,  the  order  parameter  may 
not  maintain  equilibrium  with  the  instantaneous  temperature  as  defined,  for  example,  by 
the  mean  kinetic  energy  of  the  molecular  centers  of  mass.  In  such  a  case,  the  observed 
nonlinear  refraction  will  be  reduced,  to  the  extent  that  the  order  parameter,  and  thus  the 
index  of  refraction,  is  unable  to  equilibrate  during  the  pulse. 

The  free  relaxation  of  the  orientational  order  parameter  has  been  observed  directly 
[6,  21,  9]  in  the  isotropic  phase,  and  has  been  found  to  show  a  substantial  pre-transitional 
slowing  down  as  the  transition  is  approached  from  above.  This  is  in  accordance  with 
a  theory  due  to  de  Gennes.[12]  The  time  scales  involved  are  found  to  be  on  the  order 
of  hundreds  of  nanoseconds  for  temperatures  close  to  the  transition.  On  the  basis  of 
the  theory  of  de  Gennes  [12],  it  is  expected  that  a  similar  phenomenon  will  occur  in 
the  nematic  phase,  with  a  similar  slowing  down  as  the  transition  is  approached  from 
below.  Thus,  near  the  nematic-isotropic  transition,  the  effective  value  of  n2  measured  in 
the  nanosecond  regime  will  be  determined  by  a  competition  between  a  large  coefficient, 
dn/dT,  and  slow  relaxation  of  the  orientational  order,  leading  to  a  long  rise  time  for  the 
thermal  nonlinear  refraction.  In  fact,  this  effect  has  already  been  seen  experimentally  by 
Armitage  and  Delwart,[8]  who  found  that  the  rise  time  for  a  thermal  diffraction  grating 
showed  a  large  increase  as  the  nematic-isotropic  transition  was  approached  from  below. 
They  also  attributed  the  effect  to  slow  relaxation  of  S.  Therefore,  in  order  to  predict  the 
effective  value  of  n2  which  will  be  observed  with  laser  pulses  in  the  nanosecond  regime, 
and  especially  to  account  for  the  temperature  dependence  of  n2  near  the  phase  transition, 
it  will  be  necessary  to  include  the  appropriate  relaxation  dynamics  in  the  description. 

Here  we  will  show  our  results  which  illustrate  how  the  relaxation  dynamics  of  the 
orientational  order  can  be  accounted  for  in  describing  the  thermal  nonlinear  refraction  of 
a  nematic  liquid  crystal.  In  particular,  emphasis  is  placed  on  the  role  of  slowing  down 
of  the  orientational  relaxation  near  the  nematic-isotropic  transition  in  determining  the 
observed  nonlinear  response. 
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In  Sec.  7.2.2,  some  aspects  of  the  theory  of  de  Gennes  [12]  which  underlie  the  subse¬ 
quent  treatment  of  both  equilibrium  and  relaxation  phenomena  in  nematic  liquid  crystals 
are  reviewed.  In  Sec.  7.2.3  a  simple  model  for  the  response  of  the  orientational  order  of 
a  liquid  crystal  subject  to  heating  by  an  optical  pulse  with  a  Gaussian  temporal  profile 
is  presented.  This  model  is  based  on  treating  the  non-orientational  degrees  of  freedom  as 
constituting  a  “bath”  which  is  assumed  to  equilibrate  quickly  to  a  well-defined  temper¬ 
ature.  Energy  taken  up  by  the  system  from  the  optical  pulse  can  exchange  between  the 
bath  and  the  orientational  degrees  of  freedom.  The  orientational  order  is  taken  to  relax 
exponentially  toward  the  equilibrium  value  associated  with  the  instantaneous  temperature 
of  the  bath.  In  this  way,  it  is  possible  to  obtain  a  closed  system  of  equations  for  the  time 
development  of  the  order  parameter  and  the  bath  temperature.  By  defining  a  suitable 
intensity  weighted  average  of  the  change  in  index  of  refraction,  one  can  obtain  an  effective 
n2  which  is  appropriate  for  the  types  of  measurements  described  in  Secs.  3.3.1  and  4.2 
by  PalfFy  [4,  5].  This  is  illustrated  numerically  using  parameters  for  the  material  5CB. 
It  is  found  that  the  effective  n2  is  strongly  peaked  just  below  the  transition  temperature 
for  pulse  durations  greater  than  ~  0.5/ts.  However,  for  pulse  durations  in  the  nanosec¬ 
ond  regime,  slowing  down  of  the  orientational  relaxation  near  the  transition  results  in  n2 
being  nearly  independent  of  temperature.  Thus,  the  experimental  observation[5]  that  the 
nonlinear  refraction  of  5CB,  as  measured  with  7  ns  pulses  is  independent  of  temperature 
is  not  inconsistent  with  the  nonlinearity  being  thermal  in  origin. 


7.2.2  Thermal  Properties  of  Nematic  Liquid  Crystals 


According  to  the  theory  due  to  de  Gennes[I2],  the  nematic  phase  is  characterized  by  a 
symmetric,  traceless  tensor  order  parameter,  Qij  which  is  proportional  to  the  anisotropic 
part  of  the  dielectric  tensor,  In  the  cases  of  interest  to  us,  the  nematic  order  parameter 
tensor  is  found  to  be  uniaxial  in  equilibrium,  in  which  case,  Q  can  be  put  in  the  form, 


/  -5/2  0  0  \ 

Q  =  I  0  -S/2  0  (4) 

\o  os) 

where  5  is  the  scalar  order  parameter.  The  temperature  variation  of  5  can  be  obtained 
from  the  free  energy  density, 

^(5,  T)  =  ^(T)52  -  \BS3  +  \CSA  (5) 

4  U  4 


by  applying  the  condition 


(6) 
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where  the  temperature  dependence  of  the  coefficient,  A,  is  given  by 


A(T)  =  o(T  -  Tm) 


(7) 


The  free  energy,  F,  is  that  part  of  the  total  free  energy  which  is  singular  at  the  nematic 
isotropic  transition.  T  is  associated  with  the  orientational  degrees  of  freedom  of  the  long 
molecular  axes,  and  in  the  subsequent  development  it  is  assumed  that  the  other  degrees 
of  freedom  can  be  treated  as  having  a  constant  specific  heat  over  the  temperature  range 
of  interest.  The  truncation  of  Eq.  5  at  4‘*  order  is  somewhat  arbitrary,  but  it  has  become 
conventional.  This  level  of  truncation  provides  a  minimal  description  of  the  qualitative 
features  of  the  nematic-isotropic  transition.  Applying  the  condition  of  Eq.  6  to  Eq. 
reffreee  gives  rise  to  two  solution  branches  which  satisfy  the  local  stability  condition, 


a 

as 2 


>  0. 


(8) 


For  T  >  T there  is  an  isotropic  branch  with  5  =  0.  A  nematic  branch  exists  for  T  <  T f, 
where 

r'  =  r-  +  ~  (9) 


with  an  order  parameter  given  by, 

S(T)  =  S'  +  (g(T*  -  T)) 


1/2 


where 


5f  = 


2  C' 


(10) 


(11) 


The  actual  nematic-isotropic  phase  transition  occurs  at  a  temperature,  Tc  between  T* 
and  Tf  given  by, 

2  B3 

r‘  =  r  +  9  ^  ‘12> 

which  is  the  temperature  at  which  T  goes  to  zero  for  the  nematic  branch. 


In  addition  to  the  order  parameter,  another  quantity  which  can  be  obtained  from  J-, 
which  will  be  needed  later,  is  the  part  of  the  specific  heat  associated  with  the  orientational 
degrees  of  freedom, 


C0 


-t*L 

d.T* 
o TS 
2 


[2(t'-d] 


1/2 


(13) 


where  the  total  derivative  notation  has  been  used  to  indicate  differentiation  including  the 
dependence  of  5  on  T. 
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The  uniaxial  order  of  the  nematic  phase  gives  rise  to  a  uniaxial  dielectric  tensor  with 
two  eigenvalues,  one  for  fields  polarised  parallel  to  the  director  and  one  for  fields  polarised 
perpendicular  to  the  director,  which  can  be  written 

<||  =  <1  + 

<x  =  <1  ~  (14) 

These  give  rise  to  extraordinary  and  ordinary  indices  of  refraction, 

n«  =  v^ii 

«o  =  v/*I  (15) 

In  all  the  numerical  calculations  that  follow,  the  following  set  of  parameters,  which  have 
been  obtained  for  5CB,[9]  are  used. 

a  =  0.13  x  10® J/Km* 

B  =  1.6x10 V/ro3 
C  =  3.9x10 *J/m3 

T  =  34°C  (16) 

Using  this  set  of  parameters  in  Eqs.  10  and  14,  and  by  determining  the  values  of  et  and 
e2  at  the  reference  point, [13] 

T  =  34.9 °C 
Tie  =  1.632 

na  =  1.532  (17) 

the  temperature  dependence  of  the  indices  of  refraction  shown  in  Fig.  28  is  obtained.  In 
these  plots,  it  has  been  assumed  that  the  only  temperature  dependence  of  e\\  and  C||  comes 
through  the  dependence  of  S  on  T.  Actually,  61  e2  will  also  depend  on  the  temperature 
through  their  dependence  on  the  molecular  density.  However,  the  moat  important  part 
of  the  temperature  dependence  of  the  indices  of  refraction,  particularly  near  the  phaae 
transition  is  due  to  changes  of  S  with  T.  Since  the  focus  of  this  work  is  mainly  on  the 
behavior  near  the  phaae  transition,  and  on  the  dynamics  of  the  orientational  order,  the 
density  will  be  treated  as  being  constant  throughout. 

Before  proceeding,  some  comment  regarding  the  accuracy  of  the  numerical  results  is  in 
order.  It  was  remarked  earlier  that  the  truncation  of  the  free  energy  expansion,  Eq.  5,  at 
fourth  order  in  5  is  somewhat  arbitrary.  This  truncation  has  the  result  that  the  calculated 
behavior  of  5  is  not  in  very  good  agreement  with  direct  observations[14]  for  temperatures 
more  than  a  few  degrees  below  the  phase  transition.  The  indices  of  refraction  have  a 
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Figure  28:  Predicted  temperature  dependence  of  indices  of  refraction 

corresponding  departure  from  measurements.[13]  The  value  of  S  calculated  from  Eq.  10 
does  not  saturate  properly  in  the  nematic  phase,  and  so  the  resulting  indices  of  refraction 
are  changing  more  rapidly  at  low  temperatures  than  the  experimentally  measured  ones. 
It  should  be  possible  to  correct  this  by  using  a  higher  order  expansion  for  Eq.  5,  but 
since  there  are  not  any  widely  agreed  upon  parameters  for  a  higher  order  expansion,  it 
seems  preferable  to  live  with  the  limitations  of  the  fourth  order  expansion.  The  *nai^ 
focus  here  is  on  the  behavior  near  the  nematic-isotropic  transition  where  Eq.  5  should  be 
adequate,  but  it  should  be  borne  in  mind  that  the  numerical  results  are  not  quantitatively 
accurate  more  than  a  few  degrees  into  the  nematic  phase.  It  is  expected  that  in  general, 
the  qualitative  conclusions  will  be  valid. 

Fig.  28  shows  that  the  indices  of  refraction  have  a  strong  temperature  dependence 
close  to  the  nematic-isotropic  transition.  According  to  Eq.  3,  this  is  expected  to  lead 
to  a  large  effective  A).  It  has  been  argued  that  this  can  be  counteracted  for  short  pulse 
durations,  by  slow  relaxation  of  the  orientational  order  dose  to  the  transition.  Following 
DeGennes[12],  the  free  relaxation  of  the  order  parameter  is  governed  by 

■*~S  <«) 

where  v  is  a  viscosity  coefficient,  which  can  be  related  to  the  Leslie  coefficient  71  .[15]  If  Eq. 


Figure  29:  Predicted  behavior  of  relaxation  time  with  temperature 
18  is  linearized  about  the  equilibrium  state,  assuming  a  fixed  temperature,  one  obtains 

»S  =  ~(S-S*)  (19) 

where  So  is  the  equilibrium  value  of  5  at  the  given  temperature.  This  leads  to  exponential 
relaxation  of  the  order  parameter  to  the  equilibrium  value,  with  a  time  constant 

which  will  be  referred  to  as  the  isothermal  relaxation  time.  This  exponential  relax* 
ation  behavior  has  been  observed  directly  in  the  isotropic  phase  the  optical  Kerr 
effect.[6,  21,  9]  In  the  experiments,  a  strong  optical  pulse  induces  alignment  of  the 
molecules  in  an  isotropic  fluid.  The  subsequent  decay  of  the  orientational  order  is  moni¬ 
tored  by  observing  the  induced  birefringence  with  a  weak  probe  beam.  It  is  expected  that 
the  same  description  will  hold  in  the  nematic  phase,  although  it  will  be  more  difficult  to 
observe  directly.  Taking  the  value  of  v  from  Coles  [9],  and  using  the  free  energy  expansion 
given  above,  one  obtains  the  behavior  shown  in  Fig.  29  for  the  temperature  dependence 
of  rP.  In  the  isotropic  phase,  r,  is  simply  u/a(T  —  T*),  as  has  been  observed  in  exper¬ 
iments.  In  the  nematic  phase,  the  temperature  dependence  of  rr  is  more  complicated, 
but  rr  tends  to  diverge  as  (?*  —  T)1^2  on  approaching  the  transition  from  below.  The 
large  value  of  rr  near  the  nematic-isotropic  transition  is  associated  with  the  fact  that  the 
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transition  temperature  is  Terj  close  to  the  critical  temperatures,  7*  and  T*.  It  has  been 
assumed  here  that  v  is  independent  of  temperature.  There  is  some  evidence  that  v  has 
some  temperature  dependence, [6]  however  this  should  be  a  weak  correction  to  the  rapid 
variation  of  r,  near  the  phase  transition. 


7.2.3  Orientational  Relaxation 


In  treating  the  time  dependence  of  the  orientational  order,  and  hence  the  refractive  in¬ 
dices,  of  a  liquid  crystal  sample  subjected  to  heating  by  an  optical  pulse,  the  simplification 
is  adopted  that  all  the  non-orientational  (vibrational,  translational,  etc.)  degrees  of  free¬ 
dom  constitute  a  “bath”  which  equilibrates  quickly  to  a  well-defined  temperature,  T. 
This  seems  to  be  a  reasonable  assumption,  since  the  time  scales  of  interest  here,  tens 
to  hundreds  of  nanoseconds,  are  long  compared  to  the  typical  time  scales  for  relaxation 
processes  in  liquids.  This  bath  is  assumed  to  be  characterised  by  a  specific  heat,  £7b,  so 
that  the  balance  of  energy  is  given  by, 

+  (21) 

where  a  is  the  linear  absorption  coefficient  and  I  is  the  intensity  of  the  optical  pulse,  is 
the  internal  energy  of  the  orientational  degrees  of  freedom.  In  addition,  the  orientational 
order  is  assumed  to  relax  toward  the  equilibrium  value  associated  with  the  instantaneous 
temperature,  according  to  Eq.  18.  It  is  assumed  that  the  system  is  initially  in  an  equi¬ 
librium  state,  characterised  by  QT IQS  —  0.  By  linearising  with  respect  to  the  quantities 
AT  and  A S,  which  are  assumed  to  be  small,  then  Eq.  18  becomes, 

where  all  coefficients,  such  as  Q^F/OS2  should  be  understood  as  being  evaluated  at  the 
initial  state. 

The  first  order  system,  Eqa.  21  and  22,  can  be  solved  straightforwardly,  giving 

AS(f)  =  f  F(t')  (23) 

where  F  is  the  fluence,  the  time  integral  of  the  intensity, 

F(t')  =  l*  dt"I{t")  (24) 

#-00 
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sad  the  time  cons  tint,  r*  is  related  to  the  isothermal  relaxatioa  time,  rr,  by 


r*  = 


Ck 

'Ci  +  C. 


(JS) 


Here  Ct  is  the  total  specific  heat,  which  is  the  sum  of  the  part  associated  with  the  ori¬ 
entational  degrees  of  freedom,  which  is  calculated  according  to  Eq.  13,  and  the  part,  Cs, 
associated  with  the  other  degrees  of  freedom,  which  is  assumed  to  be  temperature  inde¬ 
pendent.  The  value  used  for  £7*  is  1.78  x  10s  J/K  m3,  which  is  extracted  from  the  data  of 
Iannacchione  and  Finotdlo[16],  as  the  value  of  the  specific  heat  a  few  degrees  above  the 
transition  where  it  is  nearly  constant,  and  by  assuming  a  density  of  1.0  gm/cm3. 


In  experiments  with  pulsed  lasers,  the  optical  intensity  typically  has  a  Gaussian  tem¬ 
poral  profile, 

I(t)  =  /oe"(‘/r)*  (28) 

and  this  is  assumed  to  be  the  case  henceforth.  With  this  assumption,  the  fiueace  is, 

F{t)  =  Tlt£(l  +  erf<.tlT)).  (27) 

If  one  defines 

then  the  solution  for  AS  can  be  written  in  the  dimensionless  form, 


A5(t)  =  A50p(t/r;0  (29) 

where  (  =  t'/t  is  the  ratio  of  the  time  constant  to  the  pulse  duration,  and  the  function 
g  is  defined  to  be 

<7(2;  0  =  ^  J  ^  dxi  )/(  (1  +  er/(*i)) .  (30) 

It  is  straightforward  to  verify  that  as  x  -*  00, 


*(*.0  -*  1 


(31) 


so  that  as  t  -*  00,  AS(t)  — »  A50,  independent  of  the  value  of  the  time  constant.  It  is 
also  easy  to  show  that  for  (  «  1,  that  is  for  small  values  of  the  time  constant, 

+  *’/(*))  (32) 


so  that  in  this  case  A S{t)  is  simply  proportional  to  F(t),  reproducing  the  quasi-static 
limit.  Using  Eq.  21  the  time  dependence  of  the  temperature  can  also  be  expressed  as 


AT(<)  1  (F(t)  ,AS(t)\ 

AT0  r  ^  f„  11  >  AS0  ) 


(33) 
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Figure  30:  Time  Dependence  of  order  parameter  and  temperature. 


where  r  =  C^/C%  and 

Fo  =  y/raloT 

A T.  -  («) 

Eqs.  26, 27,  29  and  33  are  illustrated  numerically,  for  various  values  of  ( in  Fig.  30.  The 
value  of  r  =  Ct/Ct  is  needed  to  define  the  temperature.  The  value  r  =  0.5  is  a  typical 
value,  and  is  used  in  all  the  plots  of  Fig.  30  for  illustrative  purposes.  In  the  plots,  the 
fluence,  temperature  and  order  parameter  are  all  normalised  so  that  they  approach  1  as 


t  — »  oo.  In  Fig.  30a,  the  value  of  (  =  t J/r  is  0.1,  so  that  the  time  constant  is  shorter  than 
the  pulse  duration.  In  this  case  both  the  temperature  aid  the  order  parameter  follow  the 
fluence  very  closely,  and  the  system  is  nearly  in  equilibrium  at  all  times.  Fig.  30b  shows 
the  case  f  =  2.0  so  that  the  time  constant  is  twice  as  lot$  as  the  pulse  duration.  In  this 
case  there  is  a  very  noticeable  departure  from  equilibrium.  The  order  parameter  is  not 
able  to  keep  pace  with  the  fluence,  and  this  means  that  the  bath  initially  takes  up  more 
of  the  energy  than  would  be  the  case  in  equilibrium,  resulting  in  the  bath  overshooting  its 
final  temperature.  The  change  in  order  parameter,  and  hence  the  refractive  index,  during 
the  pulse  is  substantially  reduced  from  the  equilibrium  response.  The  order  parameter 
does  not  equilibrate  until  after  the  pulse  has  passed.  Fig.  30c  shows  the  case  (  —  5.0. 
This  is  qualitatively  similar  to  Fig.  30b  but  there  is  an  even  greater  overshoot  of  the  bath 
temperature  and  further  suppression  of  the  change  in  S  during  the  pulse. 

In  an  experiment,  such  as  those  described  in  Sec.  4,  in  which  the  self-refraction  of  a 
pulse  is  used  to  probe  the  nonlinear  refraction,  only  the  change  of  n  which  occurs  during 
the  pulse  will  effectively  contribute  to  the  measurement.  Since  in  the  s-scan  measurements 
the  signal  is  proportional  to  An,  and  the  weight  with  which  it  is  detected  is  proportional 
to  /,  it  is  reasonable  to  define  the  intensity  averaged  quantities 


-  JdU(t) 

=  £/*-«**«*> 

(35) 

and 

S  =  |S5. 

(36) 

Using  Eq.  29,  AS  can  be  put  in  the  form 

£3  =  ASoMO 

(37) 

where  the  function  h  is  defined  as 

h w = /_00<k*(z;*)e~*1 

(38) 

By  comparing  the  intensity  averaged  change  in  index  of  refraction 
from  an  instantaneous  nonlinear  refraction,  it  is  possible  to  define  an 
instantaneous  nonlinearity,  for  which 

with  that  obtained 
effective  »2.  For  an 

A»(<)  =  inj£J(t) 

(39) 

the  intensity  averaged  change  in  index  of  refraction  is 

-r—  2x/o 

An  =  —7= — n2 
y/2cn 

(40) 
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«• 

0.13  xlO*  J/K  m* 

B 

1.6  xlO*  J/m* 

C 

3.9  xlO*  J/m* 

T* 

34  *C 

V 

7.3  x  10"a  kg/m  • 

1.632 

1.532 

cv 

1.78x10*  J/Km* 

Table  6:  Phenomenological  parameters  used  in  calculations 

where  c  is  the  speed  of  light  and  n  is  the  linear  index  of  refraction.  Using  this  relation, 
the  effective  nonlinear  refractive  index  is  defined  in  terms  of  the  intensity  averaged  An  as 


_«//  _  V^cn-gr 
^  ~  2r/0  A 

(41) 

Together  with  Eqs.  37  and  36,  this  leads  to 

(42) 

Using  the  definition  28  and  the  relation  15,  one  obtains 

n5{/  =  2nSiifc(0 

(43) 

and 

«//  _  _£_«// 

*41  — 

(44) 

where 

0  Mj  dS  OCT 

“  4v®  &  C, 

(46) 

which  is  just  the  value  of  ity)  in  the  quasi-static  Emit. 


For  convenience)  the  phenomenological  parameters  used  in  the  calculations  are  sum¬ 
marised  in  Table  6.  Fig.  31  shows  the  calculated  results  for  along  with  n§  for  several 
values  of  the  pulse  duration.  The  value  of  the  linear  absorption  coefficient)  a,  is  taken 
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Figure  31:  Calculated  temperature  dependence  of  n2  for  three  pulse  durations,  r  = 
650,65,6.5ns.  Dashed  curves  are  the  same  in  each  figure  and  show  quasi-static  n£.  Solid 
curves  are  . 
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to  be  10  m~*.  The  actual  value  o l  a  hot  5CB  is  not  well  known,  so  this  value  of  alpha  is 
chosen  to  approximately  agree  with  measurements  for  related  materials. [17]  The  value  of 
a  enters  into  the  plots  as  an  overall  factor,  through  as  given  in  Eq.  45.  Otherwise, 
the  behavior  shown  in  the  plots  is  completely  determined  by  independent  measurements. 

Fig.  31a  shows  the  temperature  dependence  of  and  n§  for  polarisation  parallel 
to  and  perpendicular  to  the  director  for  a  pulse  duration  of  650  ns.  At  this  time  scale 
the  behavior  is  practically  quasi-static,  and  there  is  very  little  departure  of  from 
n°.  Fig.  31b  shows  the  same  result  for  a  pulse  duration  of  65  ns.  In  this  case,  the  slow 
relaxation  of  the  orientational  order  near  the  phase  transition  has  mostly  canceled  the 
enhancement  of  n2  due  to  the  large  value  of  dn/dT.  Fig.  31c  shows  the  same  thing 
again  for  r  =  6.5  ns,  which  is  the  time  scale  probed  in  the  experiments  of  Sec.  4.  In 
this  case,  there  are  significant  effects  of  the  finite  relaxation  time  over  the  whole  range 
of  temperatures  shown.  The  n\**  is  nearly  independent  of  temperature,  and  decreasing 
slightly,  rather  than  increasing  in  magnitude  on  approaching  the  transition  temperature. 
The  order  of  magnitude  of  the  calculated  TI3  is  consistent  with  the  measurements  of  Sec. 
4,  although,  given  the  uncertainty  of  linear  absorption,  a,  more  detailed  comparisons  with 
the  measured  magnitude  of  n3  are  not  meaningful. 


7.2.4  Analysis  of  Thermal  and  Relaxation 


Our  simple  formulation  accounts  for  the  relaxation  dynamics  of  the  orientational  order  of 
a  nematic  liquid  crystal  subjected  to  heating  by  linear  absorption  of  energy  from  an  optical 
pulse.  This  enables  the  calculation  of  the  time  dependence  of  the  order  parameter  and 
thus  the  refractive  indices.  By  defining  a  suitable  intensity  weighted  average  of  the  change 
in  index  of  refraction,  an  effective  nonlinear  refractive  index,  ,  is  extracted,  which  is 
comparable  to  what  would  be  obtained  from  a  measurement  in  a  z-scan  experiment  of  the 
type  described  in  Sec.  4  which  we  have  also  published,  [4,  5]. 

Our  main  conclusion  is  that  is  essentially  independent  of  temperature  for  pulses 
of  several  nanosecond  duration.  This  is  due  to  the  fret  that  the  relaxation  time  of  the 
orientational  order  shows  strong  enhancement  as  the  nematic-isotropic  transition  is  ap¬ 
proached  from  below,  and  for  short  pulses,  this  effect  counteracts  the  strong  increase  of 
dn/dT  near  the  transition. 

With  this  result,  it  seems  possible  to  explain  many  of  the  the  features  of  the  observed 
nonlinear  refraction  in  the  nanosecond  regime  in  terms  of  a  thermal  mechanism.  For 
example,  it  was  stated  in  Sec.  5.1  that  for  all  the  materials  studied,  n2y  <  0  and  n2x  > 
0,  as  would  be  expected  if  the  main  effect  were  reduction  of  orientational  order  due 
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to  heating.  In  addition,  experiments  of  Sec.  4.1  vith  two  pulses  separated  by  several 
nanoseconds  showed  a  clear  influence  of  the  first  pulse  on  the  second,  which  is  consistent 
with  the  accumulation  of  heat  in  the  sample.  The  same  type  of  temporal  behavior  was  seen 
in  the  numerical  simulation  of  pulse  propagation  in  Sec.  7.3  which  assumed  a  thermal 
mechanism.  The  magnitude  of  the  observed  nonlinearity  seems  to  be  consistent  with 
heating,  up  to  the  uncertainty  in  the  linear  absorption  coefficient. 


However,  a  complete  description  of  the  thermal  nonlinear  refraction  requires  some 
extensions  of  this  work.  The  thermal  variation  of  the  density  has  been  neglected,  which 
has  the  effect  that  the  calculated  ftj  is  sero  in  the  isotropic  phase.  In  the  measurements, 
is  found  to  be  smaller,  but  still  significant  in  the  isotropic  phase,  indicated  that  the 
density  variations  make  a  non-negligible  contribution.  The  present  calculations  satisfy 


_«//  _  _*«// 
nax  -  2n*ll 


(«) 


strictly,  whereas  the  experimental  measurements  satisfy  this  only  approximately.  This  is 
partly  due  to  the  density  variations,  which  will  make  a  negative  contribution  to  n3  for 
both  polarisations.  It  is  also  likely  that  anisotropy  of  the  optical  absorption  coefficient  a 
contributes  to  departures  from  Eq.  46.  Certain  simplifications  have  been  made  here  in 
order  to  focus  attention  on  a  specific  point,  and  in  order  to  avoid  the  introduction  of  a 
large  number  of  adjustable  parameters.  The  main  result  is  expected  to  be  independent  of 
these  simplifications. 


7.3  Nonlinear  Optical  Propagation  Modeling 


We  carried  out  a  theoretical  analysis  of  the  z-scan  experiments  used  to  measure  the 
nonlinear  optical  properties  of  liquid  crystals,  with  the  aim  of  clarifying  the  mechanism 
responsible  for  the  observed  nonlinearity  in  the  nanosecond  regime.  We  considered  various 
orders  of  nonlinear  absorption  and  nonlinear  refraction  in  order  to  examine  iheir  effects  on 
the  resulting  z-scan  curve.  The  shape  of  the  z-scan  curve  in  principle  contains  information 
about  the  order  of  the  nonlinearity.  We  found  that  qualitative  features  of  the  experimental 
z-scan  curves  are  consistent  with  a  fifth  order  absorption  and  a  third  order  refraction. 
However,  given  the  level  of  disagreement  between  the  calculated  and  experimental  curves, 
it  is  not  possible  to  make  quantitative  statements  about  the  coefficients  of  the  various 
orders.  The  magnitude  of  the  refractive  part  of  the  nonlinearity  is  consistent  with  laser 
heating  due  to  linear  absorption.  The  temporal  profile  of  a  two  pulse  experiment  is  also 
consistent  with  this  mechanism.  We  suggest  that  the  temperature  dependence  of  the 
nonlinear  refraction  might  be  explained  in  terms  of  a  thermal  mechanism  by  accounting 
for  the  temperature  dependence  of  the  rate  of  relaxation  of  the  orientational  order. 
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7.3.1  Experiments  Being  Modeled 


The  experimental  effort  provides  data  on  nonlinear  parameters  for  pulse  durations  in  the 
millisecond,  nanosecond  and  picosecond  regimes.  We  have  also  experimentally  studied 
the  dependence  of  these  nonlinear  susceptibilities  on  temperature  and  on  the  intensity  of 
the  probing  beam.  We  carried  out  these  experiments  in  geometries  in  which  reorientation 
of  the  macroscopic  director  field  is  not  expected  to  occur.  In  this  section  we  present  a 
theoretical  analysis  of  the  experiments,  focusing  attention  on  measurements  of  5CB  (4- 
cyano-4’-n-pentylbiphenyl)  in  the  nanosecond  regime  where  the  mechanism  underlying 
the  observed  nonlinear  refraction  is  not  known. 

The  s-scan  technique  is  designed  to  measure  the  nonlinear  refractive  index,  n2,  which 
is  proportional  to  the  third  order  polarizability,  *3.  However,  in  many  cases,  due  to  the 
temporally  and  spatially  nonlocal  character  of  the  response  function,  X3,  the  measured 
value  of  n2  depends  strongly  o*  the  geometrical  and  temporal  properties  of  the  measuring 
laser  beam  as  well  as  intrinsic  material  properties.  For  example,  if  the  measured  response 
is  due  to  laser  heating  of  the  sample,  we  expect 

n2  ac  An// 

dt;  nr 

~  dT  C 

where  a  is  the  linear  absorption  coefficient,  C  is  the  specific  heat  per  unit  volume,  I  is 
the  beam  intensity  and  r  is  the  pulse  duration.  In  Eq.47,  we  have  assumed  that  the  main 
source  of  heat  is  linear  absorption  of  energy  from  the  laser  beam,  and  that  the  change  in 
index  of  refraction,  n,  follows  the  equilibrium  rate  of  change  with  temperature,  dn/dT. 
Eq.  47  also  assumes  that  the  pulse  duration,  r,  is  short  compared  with  the  time  scale  for 
diffusion  of  heat  across  the  beam  profile,  which  is  given  by 

td  =  wl/D  (48) 

where  ti/0  is  the  beam  waist,  and  D  is  the  thermal  diffusion  constant.  For  values  of  D 
typical  of  liquid  crystals,  and  for  a  beam  waist  of  the  order  10  fim,  one  obtains,  77?  ~  1 
ms,  which  is  very  long  on  the  scale  of  the  measurements  we  consider  here,  for  which  r  =  7 
ns. 


(47) 


The  present  analysis  is  concerned  with  addressing  the  mechanism  behind  the  observed 
nonlinear  index  of  refraction  for  pulses  on  the  nanosecond  time  scale.  In  particular,  we 
focus  on  laser  heating  as  a  possible  mechanism.  Our  consideration  of  laser  heating  is 
motivated  by  certain  experimental  facts.  First,  the  signs  of  the  observed  ns’s  are  con¬ 
sistent  with  what  is  expected  from  laser  heating.  In  particular,  it  has  been  observed[5] 
that  in  a  large  number  of  samples,  n2  is  negative  for  polarization  of  the  beam  parallel  to 
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the  nematic  director,  while  n2  is  positive  for  polarization  per  «mdicular  to  the  director. 
This  is  what  one  would  expect  from  a  picture  in  which  the  main  effect  is  due  to  a  re¬ 
duction  of  nematic  ordering,  and  thus  optical  anisotropy,  by  heating  of  the  sample.  We 
also  note  that  other  potential  mechanisms,  namely  field  induced  orientational  order  and 
electrostriction,  predict  a  positive  n2  regardless  of  polarization.  Non-resonant  electronic 
hyperpolarizability  effects  tend  to  produce  a  positive  fi2  for  polarization  parallel  to  the 
director  and  variable  sign  for  the  perpendicular  polarization.  [20] 

Second,  the  size  of  the  observed  nonlinear  refraction  is  consistent  with  what  would  be 
predicted  given  reasonable  values  of  material  parameters.  This  is  shown  more  clearly  in 
the  following  section  where  we  show  comparisons  of  calculated  and  experimental  z-scan 
curves.  The  reported  nonlinear  refractive  index  (Sec.  4)  of  —54  x  10~u  esu  for 

5CB  measured  with  a  7  ns  pulse  is  very  large,  and  is  more  than  an  order  of  magnitude 
larger  than  the  value  obtained  with  a  33  ps  pulse.  Third,  we  note  that  Eq.  47  implies 
that  if  laser  heating  is  the  important  mechanism,  then  the  observed  n2  should  scale  in 
proportion  to  the  pulse  duration. 

In  Secs.  4  and  3  we  described  the  measurement  of  n2  for  several  nematic  liquid  crystals 
with  pulses  of  both  10  ms  and  7  ns  of  duration.  From  our  experimental  work,  it  was 
concluded  that  the  observed  n2  in  the  10  ms  case  was  due  to  laser  heating.  If  Eq.  47 
holds,  then  the  observed  n2  in  the  7  ns  case  should  be  ~  10'  times  smaller  than  that  in 
the  10  ms  case.  In  fact,  the  observed  n2  in  the  7  ns  case  for  5CB  and  for  ZLI-2303  are 
of  order  10  times  larger  than  this  prediction.  We  note,  however,  that  10  ms  pulses  are 
too  long  for  Eq.  47  to  hold.  The  response  over  10  ms  is  expected  to  be  reduced  from  the 
prediction  of  Eq.  47  by  a  factor  ~  td/t  ~  1/10,  and  it  is  possible  that  the  discrepancy 
could  be  explained  by  accounting  for  diffusion  of  heat  in  the  10  ms  case. 

Given  these  indications  in  favor  of  a  thermal  interpretation,  we  must  note  that  in 
the  case  of  the  nanosecond  pulses  there  is  a  shortcoming  of  this  mechanism,  which  is 
that  it  apparently  does  not  give  the  experimentally  observed  temperature  dependence  of 
n2.  Specifically,  the  strong  temperature  dependence  of  dn/dT  near  the  nematic-isotropic 
transition  is  expected  to  lead  to  a  corresponding  temperature  dependence  of  n2.  This  is 
observed  for  10  ms  pulses,  but  in  the  case  of  the  7  ns  pulses,  n2  is  found  to  be  essentially 
independent  of  temperature  right  up  to  the  phase  transition  temperature.  However,  we 
point  out  that  the  equilibrium  value  of  dn/dT  gives  an  upper  bound  on  the  variation  of  n 
with  T  measured  in  a  dynamical  experiment.  The  temperature  dependence  of  the  index 
of  refraction,  n,  is  mainly  due  to  the  temperature  dependence  of  the  order  parameter, 
5,  with  dn/dT  <x  dS/dT.  If  the  relaxation  time  for  variations  of  the  order  parameter, 
rr,  exceeds  the  duration  of  the  pulse,  then  it  is  expected[2]  that  the  observed  n2  will  be 
reduced  by  roughly  r/rr.  According  to  a  theory  of  de  Gennes,[12]  the  orientational  order 
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relaxes  according  to 


7 dS/dt  =  ~(S  -  So)  (49) 

where  7  is  a  viscosity  coefficient.  This  gives  r,  =  7 /fjf.  As  one  approaches  the  nematic* 
isotropic  transition  from  below,  approaches  aero  as  (T*  —  T)l/a,  where  Tf  is  the 
maximum  superheating  temperature,  typically  a  fraction  of  a  degree  above  the  nematic- 
isotropic  transition  temperature.  This  leads  to  very  large  values  of  the  relaxation  time 
near  the  transition.  At  the  same  time  dSjdT  ~  (Tf  -  T)~1^2,  so  that  for  pulses  that 
are  short  compared  to  rr,  the  effects  of  orientational  relaxation  may  cancel  the  effects 
of  rapid  variation  of  5  near  the  transition.  Slowing  down  of  the  orientational  relaxation 
has  been  observed  directly[6,  21]  in  the  isotropic  phase  using  the  optical  Kerr  effect,  in 
which  case  rr  ~  (T  —  T*)-1.  Relaxation  times  of  the  order  of  hundreds  of  nanoseconds 
have  been  seen  in  these  experiments.  In  the  calculations  that  follow,  we  will  consider  the 
measurements  on  5CB  at  25  °C,  which  is  10  °C  below  the  nematic-isotropic  transition,  so 
that  these  slowing  down  effects  may  not  be  important.  However,  any  attempt  to  describe 
the  temperature  dependence  of  the  thermal  nonlinearity  near  the  transition  must  include 
these  effects  if  it  is  to  be  valid  in  the  nanosecond  regime. 


7.3.2  Pulse  Propagation 


In  this  section  we  present  the  theory  behind  our  detailed  simulations  of  z-scan  experiments. 
The  z-scan  method  is  a  single  beam  technique  for  measuring  the  nonlinear  absorption 
and  nonlinear  refractive  index  of  a  sample.  [10, 11]  In  this  technique,  the  sample  is  moved 
along  the  axis  of  a  focused  beam,  and  the  pulse  energy  transmitted  through  the  sample  is 
measured  as  a  function  of  sample  position.  With  an  aperture  in  front  of  the  detector,  the 
measurement  is  sensitive  to  both  nonlinear  refraction  and  nonlinear  absorption.  Without 
an  aperture,  or  with  a  large  aperture,  the  total  transmitted  energy  is  measured,  sensing 
only  the  nonlinear  absorption.  Usually,  both  types  of  measurements  are  performed,  the 
large  aperture  scan  giving  an  independent  measurement  of  the  nonlinear  absorption.  In 
our  calculations,  we  assume  that  the  incident  pulse  is  represented  by  the  fundamental 
gaussian  beam,  [22] 

Ein{r} z. , <)  =  E0^~r  exp(-i2/2r2) exp(»(kz,  -  $)  -  r2(l/w(z,)2  -  ik/2R{z.)))  (50) 

W(Zi) 

where  w0  is  the  beam  waist,  and 

4>(zt)  =  arctan  (zt/zo)  (51) 

w2{z.)  =  w2  |l  +  J  (52) 
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where  zq,  the  diffraction  length,  if  given  by 


(M) 


<o  =  lni/2  (54) 

end  k  is  the  vacuum  wavenumber  for  the  beam.  We  use  the  symbol  i,  to  denote  the 
position  of  the  sample  along  the  beam  axis,  relative  to  the  focus. 

Given  the  field  at  the  incoming  surface  of  the  sample,  the  effect  of  propagation  through 
the  nonlinear  medium  is  to  attenuate  the  bea-  distort  the  wave  fronts  through 

nonlinear  refraction.  In  the  limit  of  a  thin  sam^  *o,  these  effects  are  represented 

by  the  equations, 


dlfdz  =  -(a  +0I  +  id7)I  (55) 

dA^/ds  =  A  An  (56) 

where  a,0  and  fi  are  the  linear,  third  order  (two  photon)  and  fifth  order  (three  photon) 
absorption  coefficients.  An  is  the  change  in  index  of  refraction  due  to  nonlinear  effects, 
and  in  the  present  calculation,  An  is  modeled  by, 

An(r,  2,  t)  =  7 1(r,z,t) 

+dn/dT~  J  ^  di!  [/ao/(r,  2,0  +  0  +  *,  t')]  (57) 

The  first  term  in  Eq.  57  represents  a  third  order  refractive  effect  which  is  instantaneous 
on  the  time  scale  of  the  pulse.  The  second  term  represents  the  effect  of  a  change  in  tem¬ 
perature  of  the  sample,  due  to  the  deposition  of  energy  by  the  beam  through  the  various 
absorption  processes.  The  factors  ffi  and  fft  are  included  as  adjustable  constants, 
0  <  /  <  1,  to  account  for  the  fact  that  only  a  fraction  of  the  absorbed  energy  may  be 
converted  to  heat  in  the  duration  of  the  pulse. 

Having  solved  Eqs.  55-57,  the  field  at  the  exit  surface  of  the  sample  is  given  by 

E(r,  z ,  +  I,t)  =  Ein{r ,  2,,  t)y/l(r,z.  +  L,t)/Iin(r,  z„t)  exp(»  A#r,  2„  t))  (58) 

where  L  is  the  sample  thickness.  To  propagate  the  field  from  the  sample  to  the  detector, 
we  make  use  of  Kirchoff’s  integral  in  the  Fresnel  diffraction  limit.[23]  In  our  simulated  s- 
scans  we  vary  the  sample  position,  z,,  and  calculate  the  pulse  energy  through  the  detector 
by  integrating  the  intensity  at  the  detector  over  the  detector  aperture  and  over  time.  We 
note  that  the  theory  presented  here  is  not  essentially  different  from  that  elsewhere[10, 11], 
except  that  we  do  not  make  use  of  the  gaussian  decompoeition[24]  to  propagate  the 
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field  from  the  sample  to  the  detector.  We  chose  the  present  method  over  the  gauaziaa 
decomposition,  since  the  gaussian  decomposition  makes  use  of  an  expansion  of  exp(s'A^) 
in  powers  of  A^,  which  in  our  case  is  of  order  unity.  Our  method  is  convenient  in  practice, 
since  the  propagation  from  sample  to  detector  is  independent  of  any  assumption  about 
the  propagation  through  the  sample,  enabling  a  great  deal  of  flexibility  with  regard  to  the 
modeling  of  nonlinear  effects. 


7.3.3  Calculated  Results 

All  of  the  subsequent  results  were  calculated  with  the  following  parameters  for  the  laser 
beam,  chosen  to  correspond  to  the  experiments  which  we  also  performed. 


Pulse  Energy: 

Pulse  Duration  (fwhm): 
Beam  Waist: 

Vacuum  Wavelength: 


24  /iJ 

6.5  ns 

7.5  fim 
532  nm  . 


We  assume  a  sample  thickness  of  25  /im  throughout,  and  all  of  the  following  results  are 
calculated  for  polarization  of  the  beam  parallel  to  the  director. 


Figs.  32  and  33  show  the  results  of  a  simulated  z-scan.  The  small  aperture  scan  was 
calculated  with  A  =  0.01,  meaning  that  in  the  absence  of  nonlinear  absorption,  1%  of  the 
incident  light  passes  through  the  aperture  into  the  detector.  In  the  large  aperture  case,  we 
integrate  all  of  the  intensity  at  the  exit  surface  of  the  sample,  so  that  the  large  aperture 
scan  is  sensitive  only  to  nonlinear  absorption.  The  calculations  are  compared  to  the 
experimental  results  for  the  material  5CB,  which  shows  some  of  the  largest  nonlinearities 
of  the  materials  studied.  The  curves  in  Figs.  32  and  33  were  calculated  with  the  following 
material  parameters,  thought  to  be  typical[2]  for  the  materials  under  study: 


dn/dT 

~ 

-3  x  nr3  jr 1 

C 

= 

2  x  106  J/m3  K 

a 

= 

35.0  m-1 

fa 

= 

1.0 

h 

= 

0.0 

U 

= 

0.0 

7 

= 

0.0. 

(59) 
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Figure  32:  Small  Aperture  Z-scan. 


Using  fa  =  1.0,/)  =  fn  —  0,  means  that  all  of  the  energy  absorbed  through  linear 
absorption  is  converted  into  heat,  bat  none  of  the  higher  order  absorption  processes  result 
in  heating  of  the  sample  within  the  duration  of  the  pulse.  This  results  in  a  third  order 
refractive  process,  that  is,  An  a  /.  The  linear  absorption,  a,  was  adjusted  to  match 
the  peak  to  peak  variation  of  the  experimental  s-scan.  Taking  7  =  0.0  that,  in 

the  calculations,  all  of  the  nonlinear  refraction  comes  from  heating.  This  is  done  for 
illustrative  purposes  here,  in  order  to  show  that  the  resultant  magnitude  of  the  nonlinear 
refraction  is  in  reasonable  agreement  with  the  measurements.  However,  if  we  take  for  7 
the  value  measured  with  33  ps  pulses,  where  instantaneous  effects  presumably  dominate, 
the  resulting  contribution  to  the  phase  shift,  Eq.  56,  is  negligible. 

The  solid  curves  in  Figs.  32  and  33  are  for  a  purely  fifth  order  absorptive  process, 
where  we  put  in  Eq.  55  0  =  0,  p  =  3.5  x  10-M  m3/W*.  Our  consideration  of  fifth 
order  absorption  here  is  motivated  by  the  experimental  observation  that  the  /?  obtained 
by  fitting  the  data  assuming  a  third  order  absorption  in  fact  has  a  substantial  linear 
dependence  on  the  incident  intensity.  The  dashed  curves  in  Figs.  32  and  33  correspond  to 
a  purely  third  order  absorption  with  p  =  0, 0  =  2.8  x  10~9  m/W.  These  values  of  0  and 
H  were  chosen  to  match  the  minimum  of  the  large  aperture  scan.  Some  understanding  of 
the  relation  between  the  shape  of  the  calculated  curve  and  the  order  of  the  underlying 
effect  can  be  gained  by  recalling  that  one  of  the  results  in  the  literature[10, 11]  is  that  one 
can  approximately  think  of  nonlinear  absorption  and  nonlinear  refraction  as  operating 
independently.  Then  the  small  aperture  scan  curve  can  be  thought  of  as  the  product  of  a 
refractive  part,  which  has  the  peak-valley  behavior,  with  the  nonlinear  absorption,  which 
has  a  minimum  at  zt  =  0.  In  general,  a  higher  order  effect  tends  to  bring  these  features 
of  the  curve  closer  to  z,  =  0.  For  example,  it  was  found  by  the  CREOL  group[ll]  that 
for  a  third  order  refraction,  the  maximum  and  minimum  of  the  small  aperture  scan  are 
separated  by  1.7zo,  while  for  a  fifth  order  refraction,  the  separation  is  1.2zo-  The  main 
effect  of  the  fifth  order  absorption  is  seen  in  Fig.  33,  which  shows  that  the  minimum  of  the 
large  aperture  scan  is  more  concentrated  near  the  origin,  in  contrast  with  the  third  order 
absorption.  In  the  small  aperture  scan,  Fig.  32,  the  refractive  mechanisms  for  the  two 
curves  are  the  same.  However,  for  the  dashed  curve,  with  the  third  order  absorption,  the 
peak  at  negative  z,  which  is  normally  associated  with  a  defocusing  medium[10, 11]  (n3  <  0) 
has  been  nearly  swallowed  by  the  broad  absorption  minimum.  For  the  solid  curve,  which 
is  calculated  with  a  fifth  order  absorption,  the  narrower  absorption  minimum  allows  for  a 
peak  in  the  small  aperture  curve  at  negative  We  note  that  the  experimental  data  do 
seem  to  show  such  a  peak.  We  find  that  this  feature  is  impossible  to  obtain  with  a  purely 
third  order  absorption.  The  experimental  large  aperture  scan  also  shows  some  sign  of 
the  narrower  minimum  characteristic  of  fifth  order  absorption,  although  this  is  difficult  to 
distinguish  from  the  scatter  in  the  data.  It  is  likely  that  more  than  one  order  of  nonlinear 
absorption  is  in  effect.  However,  mixtures  of  third  and  fifth  order  absorption  constrained 
to  give  the  same  minimum  of  the  large  aperture  scan  simply  interpolate  between  the 
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curve*  shown.  Since  this  does  not  seem  to  result  in  say  Vetter  fit  to  the  data,  particularly 
for  the  small  aperture  case,  it  does  not  seem  to  be  possible  to  obtain  definitive  values  for 
the  coefficients  by  fitting  a  single  s-scan  curve. 

Nonlinear  absorption  may  contribute  to  heating  of  the  sample.  The  experiments  show 
that  40%  of  the  incident  energy  is  removed  from  the  beam  by  these  processes  in  the 
situation  we  consider.  If  all  of  this  energy  were  converted  to  heat,  the  sample  temperature 
would  rise  by 

AT  «  fiIpr/C 

w  4000 Jf,  (60) 

where  Jp  is  the  peak,  on-axis  intensity.  Clearly  this  does  not  occur.  It  may  be  that  the 
energy  absorbed  by  these  higher  order  processes  remains  in  long-lived  electronic  excita¬ 
tions  for  the  duration  of  the  pulse,  or  is  carried  away  by  fluorescence.  This  is  not  known. 
If  some  fraction  of  the  energy  absorbed  by  higher  order  processes  is  turned  into  heat 
during  the  pulse,  then  fa  and  should  not  be  strictly  sero.  In  Fig.  34  we  illustrate  this 
possibility  by  calculating  a  small  aperture  scan  where  the  refractive  mechanism  is  entirely 
heating  due  to  fifth  order  absorption.  This  results  in  a  seventh  order  refractive  effect.  We 
have  taken  fa  =  1.6  x  10~3.  The  result  of  this  calculation  is  shown  as  the  dashed  line  in 
Fig.  34  along  with  the  solid  curve  of  Fig.  32  and  the  experimental  points. [19]  The  higher 
order  refraction  results  in  the  peak  of  the  curve  being  pulled  in  closer  to  x,  =  0,  where  it 
is  almost  completely  eliminated  by  absorption.  This  leads  to  worse  agreement  with  exper¬ 
iment,  particularly  in  the  region  of  negative  Since  it  was  reported  experimentally  [5] 
that  the  measured  n?  is  independent  of  the  input  intensity,  it  appears  that  the  dominant 
refractive  effect  is  truly  second  order. 

With  regard  to  the  thermal  mechanism,  we  note  that  although  we  have  adjusted  the 
value  of  a  used  in  the  calculations  of  Figs.  32  and  33  to  fit  the  experimental  data, 
the  resultant  value  of  a  =  35.0  m-1  is  well  within  what  is  considered  reasonable,  the 
precise  value  of  a  being  difficult  to  measure.  [2]  From  the  comparisons  of  our  calculations 
with  the  experimental  data,  we  conclude  that  the  magnitude  of  the  observed  nonlinear 
refraction,  as  measured  by  the  peak-peak  variation  of  the  s-scan  curve,  is  consistent 
with  this  mechanism.  However,  the  main  feature  of  the  thermal  mechanism  is  the  way 
the  nonlinear  refraction  accumulates  in  time,  as  in  Eq.  57.  Since  the  s-scan  technique 
measures  the  integrated  pulse  energy,  no  information  is  obtained  about  the  temporal 
character  of  the  mechanism.  Our  s-scan  curves  could  have  been  calculated  by  putting 
fa  =  fa  =  fa  =  0  in  Eq.  57  and  using  7  =  -2.7  x  10“ 18  m2/W  with  nearly  identical 
results.  In  order  to  gain  information  about  the  temporal  behavior  of  the  nonlinearities,  we 
have  performed  experiments  (Sec.  4.1  and  Ref.  [5]  )  in  which  we  observed  the  temporal 
profile  of  the  outgoing  pulse  with  two  closely  spaced  input  pulses.  Fig.  35  shows  a 
simulation  of  just  such  a  situation. 
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Figure  34:  Z-scan  resulting  from  purely  5th  order  heating. 
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Figure  35:  Two  pulse  Z-scan. 

Temporal  profile  of  a  double  pulse.  The  solid  curre  is  for  a  sample  position  far 
from  the  focus.  The  dashed  curve  corresponds  to  a  sample  position  near  the  maximum  of 
the  small  aperture  s-tcan 


The  material  parameters  are  the  same  as  those  used  for  the  solid  curve  of  Fig.  32, 
with  the  nonlinear  refraction  again  being  assumed  to  be  entirely  thermal.  The  solid  curve 
of  Fig.  35  shows  the  outgoing  power  through  the  small  aperture  for  a  case  in  which  the 
sample  is  located  far  from  the  beam  focus,  so  that  nonlinear  effects  are  not  important.  In 
this  case  the  output  power  is  simply  proportional  to  the  input  power.  The  dashed  curve 
shows  the  case  in  which  the  sample  is  positioned  near  the  maximum  of  the  small  aperture 
scan,  Fig.  32.  In  this  case  there  is  important  nonlinear  refraction,  which  tends  to  increase 
the  energy  flux  through  the  aperture.  This  nonlinear  refractive  effect  builds  up  in  time 
according  to  Eq.  57  with  the  accumulation  of  heat  in  the  sample.  Since  there  is  effectively 
no  diffusion  of  heat  over  the  time  scales  we  are  considering,  the  heat  accumulated  during 
the  first  pulse  is  still  in  place  during  the  second  pulse,  leading  to  a  strong  influence  of  the 
first  pulse  on  the  second.  This  type  of  behavior  has  been  observed  experimentally.  [5] 


7.3.4  Interpretations 

We  have  performed  numerical  simulations  of  z-scan  experiments  for  comparison  with  the 
experimental  results  published  [4,  5,  19]  and  described  in  Sec.  4.  In  particular,  we  have 
considered  the  z-scan  experiments  with  7  ns  pulse  duration  and  the  material  5CB,  which 
was  the  most  extensively  studied.  We  have  shown  calculations  to  illustrate  the  z-scan 
signature  associated  with  various  orders  on  nonlinear  absorption  and  nonlinear  refraction. 
Qualitative  features  of  the  experimental  z-scan  curves  seem  to  be  consistent  with  third 
order  refraction  and  fifth  order  absorption.  We  find  that  the  magnitude  of  the  observed 
nonlinear  refraction  is  consistent  with  laser  heating  of  the  sample  due  to  linear  absorption. 
We  have  also  shown  a  calculation  of  the  temporal  profile  of  the  detected  power  through  a 
small  aperture  in  the  case  where  the  incident  beam  consists  of  two  closely  spaced  gaussian 
pulses.  Due  to  the  cumulative  nature  of  the  thermal  mechanism,  the  presence  of  the  first 
pulse  has  a  strong  influence  on  the  second.  Our  calculated  result  is  qualitatively  similar 
to  the  our  experimental  measurements,  which  clearly  shows  such  a  cumulative  effect. 

In  our  calculations,  we  have  only  considered  one  order  of  nonlinear  refraction  or  non¬ 
linear  absorption  at  a  time.  In  reality,  various  orders  of  nonlinearity  will  be  operative 
simultaneously.  In  principle,  the  shape  of  the  z-scan  curve  contains  information  about 
the  coefficients  of  various  orders.  In  practice  we  find  the  level  of  disagreement  between 
calculated  and  experimental  curves  precludes  any  such  detailed  analysis.  In  addition  to 
random  scatter  of  the  experimental  points,  which  presumably  results  from  fluctuations 
in  laser  pulse  energy,  there  seem  to  be  systematic  deviations  between  calculation  and 
measurement.  One  possible  source  of  such  deviations  is  imprecise  location  of  the  origin, 
z,  =  0,  in  the  experiments.  We  have  not  attempted  to  correct  for  this  in  the  plots  shown 
here.  Instead  of  trying  to  obtain  a  very  refined  curve  shape  for  a  single  pulse  energy, 
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it  would  seem  to  be  fruitful  to  try  to  obtain  the  coefficients  of  the  various  processes  by 
fitting  s-scan  curves  measured  over  a  range  of  pulse  energies. 

Energy  dependent  measurements  were  experimentally  measured  and  the  data  were 
analysed  to  obtain  pulse  energy  dependent  “effective”  values  for  n3  and  0.  The  effective 
T»j  was  found  to  be  independent  of  intensity,  while  0  was  found  to  depend  linearly  on 
intensity,  indicating  the  presence  of  fifth  order  absorption.  A  rough  estimate  of  the  fifth 
order  absorption  coefficient  can  be  obtained  from  their  plot  using  fi  ~  A/9/AJ  w  10~21 
m3/W2.  This  compares  favorably  with  the  value  3.5  x  10-M  m3/WJ  used  in  the  plot  of 
Fig.  32.  More  detailed  fitting  of  energy  dependent  z-scans  to  obtain  definitive  values  of 
the  coefficients  would  seem  to  be  very  profitable. 

Both  the  magnitude  and  the  temporal  behavior  of  the  observed  nonlinear  refraction 
are  consistent  with  a  thermal  mechanism.  We  have  noted  that  the  thermal  mechanism 
has  one  shortcoming  with  respect  to  explaining  the  observed  n?,  and  that  is  that  the  n2 
observed  with  7  ns  pulses  is  found  to  be  nearly  independent  of  temperature,  [5]  while  on 
the  basis  of  Eq.  47  one  expects  that  n2  will  have  the  temperature  dependence  of  dn/dT, 
which  has  rapid  variation  near  the  nematic-isotropic  transition.  However,  we  point  out 
that  7  ns  is  actually  much  shorter  than  the  characteristic  relaxation  time,  rr,  associated 
with  variations  in  the  orientational  order  parameter,  S ,  for  temperatures  near  the  tran¬ 
sition.  Given  that  rr  has  the  same  dependence  on  the  temperature  as  dS/dT  near  the 
phase  transition,  it  is  quite  likely  that  the  observed  lack  of  temperature  dependence  of 
n2  is  a  result  of  the  finite  response  time  of  5.  We  are  currently  investigating  molecular 
orientational  dynamics  theoretically,  and  should  soon  be  able  to  make  more  definite  state¬ 
ments  regarding  the  dynamics  of  the  response  of  the  orientational  order  to  variations  of 
temperature  and  external  fields,  and  the  resulting  temperature  dependence  of  n2. 

Further  experimental  data  would  also  be  very  valuable  in  identifying  the  underlying 
mechanism  of  the  nonlinear  optical  response.  For  example,  we  have  noted  that  the  cu¬ 
mulative  temporal  behavior  of  the  nonlinear  refraction  is  characteristic  of  the  thermal 
mechanism.  However,  one  can  imagine  other  mechanisms  which  would  exhibit  similar  be¬ 
havior.  Excited  state  refraction  in  which  the  molecular  polarizability  is  changed  through 
a  one  photon  excitation  would  have  a  very  similar  signature.  In  fact,  if  the  relaxation 
time  of  the  excited  state  were  long  compared  with  the  pulse  duration,  the  theoretical 
description  would  be  essentially  identical  with  Eq.  57  except  that  the  parameters  would 
be  interpreted  differently.  In  the  case  of  the  thermal  mechanism,  we  know  that  the  decay 
time,  which  is  essentially  the  diffusion  time,  Eq.  48,  is  indeed  long.  Therefore,  if  the 
thermal  mechanism  is  the  operative  effect,  one  should  not  observe  any  diminution  of  the 
effect  of  the  first  pulse  on  the  second  in  a  double  pulse  experiment  as  the  delay  between 
pulses  is  increased  to  quite  large  values  (~  1  ms).  This  is  a  quite  definite  prediction  of 
the  thermal  mechanism  and  any  departure  from  this  behavior  would  be  strong  evidence 
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•gainst  •  thermal  interpretation.  Observations  of  fluorescence  would  help  to  determine 
the  fate  of  the  large  amounts  of  energy  removed  from  the  beam  by  nonlinear  absorption. 
Also,  the  observation  of  fluorescence  at  wavelengths  shorter  than  A/2  would  confirm  the 
presence  of  fifth  order  (three  photon)  absorption. 
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Materials  Synthesis,  Acquisition  and  Preparation 

8.1  Materials  Preparation 

The  initial  throat  of  this  contract,  aa  described  in  the  proposed  work,  waa  the  investigation 
of  novel  composite  ayatema  offering  poaaibiHtiea  for  enhanced  nonlinear  effects.  Subsequent 
changes  in  the  immediate  interests  of  the  overall  DARPA  OPL  effort  resulted  in  more 
attention  to  pnre  materials. 

The  primary  activity  in  the  materials  program  for  the  first  portion  of  the  project  waa 
the  preparation  and  characterisation  of  polymer  dispersed  liquid  crystal  (PDLC)  films. 
PDLC  films  are  made  through  one  of  a  variety  of  processes.  While  considerable  experience 
exists  in  house  in  PDLC  preparation,  the  process  changes  for  each  set  of  new  materials. 
In  principle,  any  one  of  thousands  of  liquid  crystal  materials  can  be  combined  with  a 
similarly  large  number  of  polymers  to  produce  a  PDLC  film.  Since  the  preparation  condi¬ 
tions  determine  the  properties  of  the  final  film,  there  is  considerable  experimentation  and 
characterization  necessary  in  order  to  achieve  desirable  characteristics  in  any  particular 
polymer/liquid  crystal  system. 

PDLC  were  formed  either  by  polymerisation  induced  phase  separation  (PIPS)  or 
by  a  combination  of  solvent  induced  and  thermally  induced  phase  separation  (SIPS  and 
TIPS  respectively).  The  droplet  size  was  controlled  by  either  the  curing  temperature  for 
PIPS  or  by  the  rate  of  cooling  from  the  isotropic  melt  for  SIPS/TIPS. 

The  droplet  formation  process  was  studied  for  the  polymers  PMMA,  PVP  and  the 
epoxy  matrix  made  using  Epon  282,  MK107  and  Capcuxe.  In  the  initial  studies  of  PDLCs, 
the  commercially  available  liquid  crystal  E7  was  used  as  a  standard  test  case.  For  PMMA 
and  PVP  the  droplet  size a  function  of  the  rate  of  cooling  from  the  isotropic  melt.  The 
PDLC  films  were  formed  by  solvent  casting  a  chloroform  solution  on  a  glass  substrate. 
The  resulting  filing  were  further  processed  by  bong  heated  above  the  melting  point  of  the 
polymer  resulting  in  formation  of  a  homogeneous  solution.  The  films  were  then  cooled  at 
a  prescribed  rate  and  the  resulting  droplet  size  measured  using  optical  microscopy.  The 
same  film*  could  be  thermally  cycled  several  times  to  change  the  droplet  size.  For  the 
epoxy  systems,  the  cure  rate  and  concentration  of  liquid  crystal  determined  the  droplet 
size  and  concentration.  In  all  cases,  the  thickness  of  a  film  is  determined  when  it  is  sand¬ 
wiched  (prior  to  droplet  formation)  between  two  glass  substrates.  The  glass  substrates 
are  separated  by  spaces  to  produce  films  of  thickness  typically  between  10  pm  and  100 
pm. 
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The  film  preparation  supported  both  the  nonlinear  optical  measurements  activities  and 
the  IR  Shutter  activities.  For  the  IR  shutter  effort,  the  principle  effort  was  at  producing 
and  evaluating  a  variety  of  PDLC  systems.  We  had  no  experience  with  the  IR  properties 
of  these  materials  so  a  very  broad  screen  process  was  initiated.  Further  details  on  film 
preparation  specific  to  the  IR  project  are  contained  in  Sec.  9.  The  very  particular 
requirements  related  to  response  time  and  droplet  morphology  as  well  as  optical  properties 
required  a  more  extensive  investigation  into  PDLC  formation. 

The  studies  of  nonlinear  response  of  composites  relied  largely  on  the  Epon/Capcure 
epoxy  with  the  E7  liquid  crystal  during  the  first  phase  of  the  effort.  This  is  a  system  with 
which  we  had  experience  switching  in  electro-optics  applications,  so  it  required  less  initial 
investment  in  materials  preparation  effort.  The  Active  Feedback  studies  described  in  Sec. 
3.1  used  preparation  techniques  which  had  been  tested  in  liquid  crystal  display  (LCD) 
applications.  The  variable  requiring  adjustment  was  the  opacity  on  the  scattering  state, 
since  the  electronic  circuits  could  be  adjusted  to  balance  feedback  with  transmission. 


8.2  Nonlinear  Materials  Selection 


In  the  second  year  of  the  project,  a  change  in  emphasis  occurred  from  investigating  novel 
system  to  demonstrating  large  materials  nonlinearity.  This  was  a  change  in  direction 
communicated  to  all  contractors  during  the  second  annual  review  meeting. 

To  help  separate  and  identify  the  operative  mechanisms  for  OPL,  we  systemized  our 
selection  of  the  class  of  materials  to  be  studied.  The  following  properties  and  selection 
criteria  were  applied  in  our  literature  search. 

1.  A  solid  to  nematic  or  smectic  transition  should  occur  near  or  below  room  temper¬ 
ature  and  have  a  range  of  at  least  10  degrees  C  or  more  to  facilitate  experimental 
measurement. 

2.  The  material  should  be  chemically  stable  to  eliminate  uncertainties  about  contami¬ 
nates  contributing  to  or  masking  nonlinear  effects.  Similarly,  laser  induced  reactions 
would  unnecessarily  complicate  analysis  at  this  stage. 

3.  Transparency  in  the  IR  and  visible  range  is  highly  desirable. 

4.  Chemical  breadth  was  desired.  We  sought  compounds  with: 

(a)  Varying  sizes  of  conjugated  cores 

(b)  Varying  linkages  between  benzene  rings 
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(c)  Varying  substituent  groups 

(d)  Ease  of  synthesis  or  even  possible  commercial  availability. 

(e)  Simple  and  small  conjugation  structures  to  facilitate  theoretical  calculations 


Except  for  requirements  that  a  nematic  or  smectic  exist  at  a  temperature  accessible  to  the 
laser  measurement  facility,  none  of  the  above  criteria  were  deemed  absolute.  We  compiled 
a  list  of  materials  which  satisfy  all  or  most  of  the  above.  The  list  is  provided  in  Appendix 
A. 


A  subset  of  the  list  of  potentially  interesting  liquid  crystals  was  actually  chosen  for 
acquisition  and  study.  The  liquid  crystal  materials  which  were  selected  for  study  are 
listed  below.  All  were  acquired  or  synthesized.  Those  synthesized  under  this  contract  are 
indicated  as  KSU  synthesis. 


•  6DA8CB  (KSU  Synthesis) 


R  =  -02C(CH2)8CsC— C=C— C6H13 
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•  8CB  (K24),(4-cyano-4’-n-octylbiphenyl) 


•  5CB  ( 4-  cy  ano-4  n- pentylbipheny  1) 


•  E7 

Commercial  Mixture,  Proprietary 


•  OPL-7-1  (8YCB)  (KSU  Synthesis) 

4-(l-octynyl)-4'-cyanobiphenyl 


•  OPL-10-1  (5AZXB)  (KSU  Synthesis) 

4-4'-dipentylazoxy  benzene 
0 
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•  OPIrlO-2  (6AZXB)  (KSU  Syatncai*) 


4-4'-dibexyUzoxy benzene 

O 

»UC*— 


•  ZLI  2303  (PYP-606)  Phenylpyrimidine 


•  ZLI  1538  (CCH-4)  Cyclohexylcyclohexane 


CB15 


CHi  N - 
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•  T15  ( 4-cysno-4”  -  n-  alkyl-  p-  terphenyl) 


•  AANB-5  4-(ome^a-aininoaIkyloxy)-4,-nitrobiphenyl«  (KSU  Synthesis) 

•  AANB-7  4-(omepa-aminoalkyloxy)-4’-nitrobipheny!s  (KSU  Synthesis) 

•  AANB-10  4- ( omega-amino*lkyloxy )-4 nitro biphenyl*  (KSU  Synthesis) 

^TO-OH0H0)’NOI 


4-((0-aminoalkyloxy)-4'-nitrobiphenyls  (AANB) 


X  =  5,7,10 

4,4’-l2is(acryloyloxy)j2iphenyl  (BAB) 
CH2=CHCC>2— ^)h^)-o2cch=ch2 


(KSU  Synthesis) 


(KSU  Synthesis) 

4,4,-bis[^-(acryloyloxy)hexyloxy3-l,l,-l2iphenyl  (BAB-6)  v 
CH2=CHC02(CH2)60~-^^-^^--O(CH2)602CrH=CH2 
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(KSU  Syntheab) 

4,4’-bis  { 4-[6-(acryloyloxy)hexy  loxy]benzoate }  -1 , 1  '-biphcny  lcne  (BABB-6) 


CH2^HC02(CH2)60-^^C0r^^-^^-02C--^^-~0(CH2)602CCH-CH2 

(KSU  Synthesis) 

4,4’-bis  {4-[6-(methacryloyloxy)hexyloxy]benzoate  }-l,r-biphenylene  (BMBB-6) 


CH2=CCH3C02(CH2)6C>--^^-C0f^^-^^“02C-^^~0(CH2)602CCCH3-CH2 


The  specific  measurements  and  nonlinear  parameters  (or  IR  parameters)  are  described 
elsewhere  in  the  text  of  the  pertinent  experimental  section.  A  tabulation  of  materials  and 
nonlinear  parameters  are  collected  in  Appendix  A. 
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9  IR  Shutter 


Included  in  our  original  contract  was  a  separate  project  to  evaluate  the  use  of  liquid  crystal 
composites  in  infrared  shattering  applications.  The  particular  applications  were  night 
vision  devices  which  currently  employ  a  mechanical  shutter  as  a  modulator  of  incident  IR 
radiation.  The  IR  detectors  which  form  the  image  require  modulation  of  the  IR  radiation. 
We  were  not  informed  of  the  details  of  operation  or  geometry  of  current  devices  employed 
by  DoD.  We  were  however  given  some  guidelines. 

The  device  requirements  as  understood  at  the  time  of  the  proposal  were: 

1.  Acceptance  angie  of  40  degrees  -  full  angle 

2.  Electrically  line  addressable 

3.  Line  dimensions  of  0.5  inches  long  and  100  p  wide 

4.  0.5  inches  by  0.5  inches  array  Size 

5.  Low  insertion  loss,  less  than  20% 

6.  >  95%  of  incident  radiation  forward  scattered  into  an  annulus  of  radius  equivalent 
to  10  line  widths 

7.  Submillisecond  switching  times 

8.  Broadband  response  over  the  8-14  fim  region. 

The  project  developed  through  several  phases. 

1.  Concept  testing  to  determine  if  any  shuttering  of  IR  radiation  was  possible  with 
existing  PDLCs. 

2.  Materials  testing  for  contrast  enhancement. 

3.  Prototype  shutter  construction  for  transmission  enhancement. 

4.  IR  camera  device  testing 

5.  Materials  testing  for  response  time  and  scattering  angles. 


We  will  follow  the  above  outline  as  a  final  report  on  the  IR  shutter  segment  of  our  effort 
and  end  with  an  evaluation  discussion. 
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0.1  Concept  Testing 


In  the  first  quarter  of  the  effort,  existing  materials  that  are  used  in  PDLC  display  ap¬ 
plication  were  used  to  form  large  droplet  films  for  concept  testing.  Visible  wavelength 
scattering  utilised  a  polymer  film  with  liquid  crystal  droplets  in  the  .5*1  pm.  The  rule 
of  thumb  that  effective  scattering  requires  droplets  comparable  to  the  wavelength  was 
applied  to  create  PDLC  films. 

A  PDLC  film  is  prepared  by  cooling  a  solution  of  E7  in  the  thermoplastic  PMMA 
melt.  The  sise  of  the  E7  droplets  is  governed  by  the  rate  of  cooling  and  the  relative 
concentration  of  PMMA  and  E7.  We  found  the  best  results  were  achieved  with  very  large 
droplets,  «  20 pm  in  diameter  and  with  a  ratio  of  E7:PMMA  of  4:1.  The  droplets  were 
formed  by  slowly  cooling  the  solution  from  75  degrees  C  to  21  degrees  C  in  about  12 
hours,  yielding  droplets  of  20  •  30 pm  in  diameter.  A  10pm  thick  film  of  the  material  was 
sandwiched  between  NaCl  substrates  with  a  platinum  conducting  electrode.  The  platinum 
electrodes  absorb  about  50%  of  the  incident  IR.  The  IR  transmission  of  the  shutter  was 
measured  with  a  Sargent  Welch  3  -  200  IR  spectrophotometer.  The  resulting  spectra  with 
and  without  an  applied  voltage  are  shown  in  Fig.  36 

The  transmitted  intensity  is  not  quantitative  for  these  films,  however,  without  an 
applied  voltage  the  scattering  intensity  in  the  2.5  •  5  pm  region  is  on  the  order  of  the 
strongest  absorption  of  the  film  and  is  switched  to  a  highly  transmitting  state  upon 
application  of  a  voltage.  Optoelectronic  modulation  of  the  transmitted  intensity  is  seen 
at  all  wavelengths  where  there  is  no  significant  absorption;  2.5  -  3.5pm,  4.0  -  6.0pm,  9.5 
•  10.5  pm  and  8  -  14pm. 
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Figure  36:  IR  Tranomiosion  of  PMMA/E7  PDLC 
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9.2  Contrast  Enhancement 

The  effective  modulation  of  IR  radiation  by  a  PDLC  device  it  not  possible  if  there  is 


1.  High  absorption  loss, 

2.  High  surface  reflection  loss, 

3.  Poor  matching  of  the  IR  index  of  refraction  of  the  polymer  and  the  ordinary  index 
of  refraction  of  the  liquid  crystal,  or 

4.  Poor  birefringence  or  small  differences  in  the  IR  ordinary  and  extraordinary  indices 
of  refraction  of  the  liquid  crystal. 


While  there  are  other  parameters  affecting  device  performance,  the  above  were  deemed 
important  enough  to  address  first. 

The  initial  phase  of  the  research  involved  screening  of  polymers  and  liquid  crystals  for 
large  absorption  windows  in  the  2.5  -  5  fim  or  8  -  14/im  region.  The  polymers  screened 
include:  polymethlmethacrylate,  polyvinylfluoride,  polystyrene,  polyvinylacetate,  poly¬ 
carbonate,  polyvinylpyrrolidone,  and  polyvinylchloride.  Epoxy  resins  and  cure  agents 
were  also  screened.  The  polymers  all  have  absorption  bands  in  the  8  -  14 pm  region;  how¬ 
ever,  several  of  the  polymers  have  broad  transparent  regions  so  that  total  transmission  is 
acceptable  for  further  experimental  development. 

Table  7  lists  the  polymers  and  liquid  crystals  surveyed  for  use  in  this  project.  The 
IR  transmission  spectrum  of  each  of  the  pure  materials  was  obtained  to  determine  their 
absorption  characteristics  in  the  2.5  -  5  and  8  -  14 ft.n  regions.  Materials  with  broad 
transmission  windows  in  these  regions  were  used  for  formation  of  PDLC  films. 

The  E7  liquid  crystal  mixture  is  representative  of  the  IR  absorption  characteristics 
of  a  wide  variety  of  commercially  available  eutectic  liquid  crystal  mixtures  based  on 
cyanobiphenyls  and  terphenyls.  The  CM-1644  liquid  crystal  is  a  not  conjugated  liquid 
crystal  based  on  cyclohexanes  and  has  a  broad  transmission  window.  However,  it  has  a 
very  low  dielectric  anisotropy  and  therefore  requires  unacceptably  high  driving  voltages 
to  switch  to  the  transparent  state. 

PDLC  film*  were  formed  either  by  polymerization  induced  phase  separation  (PIPS)  or 
by  a  combination  of  solvent  induced  and  thermally  induced  phase  separation  (SIPS  and 
TIPS  respectively).  The  droplet  size  was  controlled  by  either  the  curing  temperature  for 
PIPS  or  by  the  rate  of  cooling  from  the  isotropic  melt  for  SIPS/TIPS.  The  IR  transmission 
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Polymer  /  liquid  crystal 

IR  Transmission 

polymethylmethacrylate  (PMMA) 

acceptable 

polyvinylpyrrolidone  (PVP) 

acceptable 

polystyrene 

acceptable 

polyvinylfluoride 

unacceptable 

polyvinylacetate 

unacceptable 

polycarbonate 

unacceptable 

polyvinylchloride 

unacceptable 

epoxies  (Epon  828,Mkl07,WC68 
with  Capcure  3*800) 

unacceptable 

E7 

acceptable 

CM*  1644 

acceptable 

Table  7:  IR  Transmission  characteristics  of  selected  polymer  and  liquid  crystals. 
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Cooling  Rate 
(°C/min) 

droplet  Diameter 
(pm) 

E7/PVP 

E7/PMMA 

>>10 

<<5 

<  <5 

4 

10-15 

3-5 

2 

15-25 

5-10 

1 

20-30 

10-15 

in 

o 

25-40 

15-25 

Droplet  size  as  a  fucntion  of  cooling  rate  for  E  7/P  VP  and 
E7/PMMA  PDLC  films. 


Table  8:  Droplet  size  as  a  function  of  cooling  rate  for  E7/PVP  and  E7/PMMA  PDLC 
films. 

spectra  of  the  resulting  PDLC  films  were  measured.  The  E7/PVP  and  the  E7/PMMA 
PDLC  films  had  the  highest  contrast  and  broadest  transmission  windows  in  the  ON  state. 

The  droplet  formation  process  was  studied  for  E7/PMMA  and  E7/PVP  PDLC  systems. 
The  droplet  size  is  a  function  of  the  rate  of  cooling  from  the  isotropic  melt.  The  PDLC 
films  were  formed  by  solvent  casting  a  chloroform  solution  on  transparent  conducting 
electrodes  with  26  or  10  pm  spacers  to  control  thickness.  The  resulting  films  were  heated  to 
above  the  melting  point  of  the  polymer  resulting  in  formation  of  a  homogeneous  solution. 
The  films  were  then  cooled  at  a  prescribed  rate  and  the  resulting  droplet  size  measured 
using  optical  microscopy.  Table  8  lists  the  cooling  rate  and  droplet  size  for  10  pm  thick 
films  of  E7/PMMA  and  E7/PVP.  The  same  films  can  be  thermally  cycled  several  times 
to  change  the  droplet  size.  The  IR  transmission  in  the  ON  and  OFF  state  was  measured 
as  a  function  of  droplet  size  for  both  PDLC  systems. 

Fig.  37  shows  the  IR  transmission  in  the  ON  and  OFF  state  of  an  E7/PVP  film 
with  a)  «  5 pm  droplets  and  with  b)  15  -  25 pm  droplets.  Comparison  of  the  OFF  state 
transmission  of  these  two  films  shows  that  smaller  droplets  results  in  higher  scattering  in 
the  2.5  •  5  pm  region  while  the  film  with  the  15  -  25  pm  droplets  has  much  higher  scattering 
in  the  9  - 14  pm  region.  This  result  agrees  with  the  general  theoretical  prediction  that  light 
scattering  will  be  most  efficient  when  the  droplet  size  is  comparable  to  the  wavelength. 
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Transalttt 


Figure  37:  Ou  and  OFF  transmission  of  small  droplet  («  5/im)  film  (top)  and  large 
droplet  (15  -  2 Opm)  film  (bottom).  Ko'"<  2  o  n+a  l  er'fii'k  aa  tn 
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Figure  38:  Experimental  Set  up  for  angular  transmission  measurements. 

9.3  Angular  Transmission 

Angular  transmission  studies  were  carried  out  using  an  experimental  set  up  shown  below. 

In  preparing  for  transmission  and  scattering  studies  of  PDLC  films  in  the  IR  when  ra¬ 
diation  is  not  at  normal  incidence,  we  measured  the  necessary  transmission  characteristics 
of  the  NaCl  coated  substrates.  These  curves  were  fit  to  determine  surface  reflectivity  and 
absorption  coefficients. 

Similar  preparatory  measurements  were  performed  on  a  mixture  of  liquid  crystal  and 
polymer  (E7  and  polyvinylpyrrolidone)  when  the  concentration  of  liquid  crystal  (E7)  was 
%  too  small  to  form  droplets.  This  material  mimics  the  polymer  host  with  its  residual 
liquid  crystal  concentration  in  the  films  where  droplets  have  been  formed.  Transmission 
studies  as  a  function  of  incident  angle  were  performed  on  these  films  without  droplets  to 
determine  their  contribution  to  background  reflection  and  absorption  profiles. 

Finally,  the  angle  dependent  transmission  characteristics  of  the  PDLC  films  were  mea¬ 
sured  at  a  number  of  points  from  2  - 18  pm  (in  regions  with  low  absorption).  The  5 pm  is 
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Figure  39:  Normalized  transmission  as  a  function  of  incident  angle  for  E7/PVP  film  with 
5,  10-15,  25-40  fim  size  droplets. 

shown  in  the  figure  as  an  illustration  of  the  results.  The  results  for  one  sample  with  80% 
E7  in  polyvinylpyrrolidone  cooled  at  different  rates  are  shown  in  the  figure  below. 

The  maximum  ON  transmission  can  be  achieved  if  there  is  perfect  index  matching 
of  the  ordinary  index  of  the  droplet,  n0,  with  the  host  polymer.  It  is  inherent  in  the 
birefringent  optical  properties  of  the  droplet  that  as  the  direction  of  the  light  becomes 
off-axis,  part  of  the  extraordinary,  n„,  index  is  encountered  during  propagation.  The 
extent  to  which  the  resulting  scattering  decreases  the  amount  transmitted  is  indicated  by 
the  difference  between  the  substrate  +  polymer  sample  and  the  substrate  +  polymer  + 
droplet  sample.  At  visible  wavelengths,  this  phenomena  manifests  itself  as  a  slight  haze 
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when  a  PDLC  is  viewed  off  axis.  Compared  to  similar  studies  in  those  systems,  the  IR 
hase  effect  is  reasonably  small. 


10  IR  Camera 

Having  produced  PDLC  films  capable  of  modulation  in  the  IR,  we  attempted  to  use  a 
liquid  crystal  composite  as  a  shutter  in  an  IR  camera.  The  first  tests  showed  only  that  it 
would  modulate  the  intensity  of  the  IR  signal  detected  by  the  camera.  Subsequently,  it 
was  shown  that  the  modulation  was  sufficient  to  produce  an  image  with  the  camera,  fig. 
40  shows  an  early  qualitative  plot  of  signal  intensity  versus  temperature  of  an  object  for 
both  mechanical  shuttering  and  PDLC  shuttering.  The  difference  between  the  brightness 
of  the  image  using  the  two  systems  shows  that  substantial  optimization  of  the  PDLC 
materials  and  the  electronics  of  the  IR  camera  remains  to  be  carried  out.  The  mechanical 
shutter  can  produce  an  image  of  an  object  which  is  only  a  few  degrees  above  background 
while  the  PDLC  begins  imaging  only  when  the  object  is  at  least  20  degrees  C  above 
background. 

We  investigated  the  effect  of  sample  thickness  on  the  contrast  of  PDLC  films  having 
relatively  large  droplets  (15  -  25  pms).  PDLC  shutters  with  a  film  thickness  of  10,  26  or 
61  pms  were  made  and  their  contrast  measured.  The  10  pm  films  show  little  modulation 
of  the  incident  radiation  whereas  the  26 pm  thick  films  switch  from  34%  to  about  52% 
transmission  between  the  OFF  and  ON  states  and  the  61/tm  thick  films  switch  between  7% 
and  21%  (note  that  the  absolute  values  for  the  percent  transmission  are  highly  dependent 
on  the  optical  arrangement  used  to  make  the  measurements,  however,  the  general  trends 
are  valid).  Much  of  the  light  lost  when  passing  through  the  PDLC  shutters  results  from 
front  and  back  surface  reflections  from  the  germanium  substrates.  We  used  antireflection 
coatings  on  the  germanium  to  increase  the  transmission  through  the  substrates  from  about 
22%  to  60%. 

We  incorporated  the  61pm  thick  PDLC  shutter  with  antireflection  coatings  in  an  in¬ 
frared  video  camera.  The  camera  with  this  shutter  could  be  used  to  image  objects  in 
our  laboratory  such  as  computer  monitors  and  human  faces.  We  did  not  measure  the 
sensitivity  of  the  camera  with  this  shutter,  however,  from  visual  observation  it  is  dear 
that  the  sensitivity  has  been  greatly  increased. 


IMAGE  BRIGHTNESS  (ARBITRARY  UNITS) 


DETECTED  IMAGE  BRIGHTNESS  AS  A  FUNCTION  OF  TARGET 
TEMPERATURE  FOR  MECHANICAL  AND  ELECTRO-OPTIC 
SHUTTERING  CONFIGURATIONS  USING  THE  MODEL  86 
THERMAL  IMAGING  SYSTEM 
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Figure  40:  Measured  image  brightness  as  a  function  of  temperature. 
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Substrate 

film 

l)  tidiness 
d  (pm) 

droplet 

diameter 

(pm) 

A  »  2-Spm 

A  a  8- 14pm  J 

%Tow 

%Torr 

C* 

*Ton 

%Topv 

C* 

(volts) 

ruff 

(ms) 

Gc 

10 

20-25 

19 

li 

23.9 

21.9 

1.1 

- 

24 

ARctd.Ce 

10 

20-25 

17.8 

1.7 

60.4 

56.1 

l.l 

29.1 

16 

ARctd.Ce 

26 

10-25 

7.4 

3.8 

51.0 

34.2 

1.5 

13.4 

254 

AR  ctd.  Ge 

61 

10-25 

3.3 

4.8 

23.0 

11.2 

2.1 

19.0 

WSi 

•  C*  ■  *Ton/*Tobi 

'Measured  for  electro-optic  response  for  A  -  8- 1 4pm 


Table  9:  Electro-optic  characteristics  of  several  PDLC  films. 

10.1  Temporal  Response 


Since  the  IR  detection  system  is  sensitive  to  changes  in  intensity  and  not  absolute  intensity, 
the  contrast  capability  of  the  IR  switch  is  important.  Furthermore,  since  there  is  a 
response  time  to  the  IR  detector,  the  temporal  modulation  characteristics  of  the  PDLC 
are  important.  Modulation  characteristics  as  a  function  of  droplet  sizes  were  investigated. 

As  expected,  the  observed  turn-off  times  for  the  PDLC  devices  are  large  (Table  9). 
Increasing  the  diameter  of  PDLC  droplets  yields  longer  turn-off  times.  The  turn-off  time 
is  defined  as  the  time  for  transmission  to  fall  from  90%  to  10%  after  removal  of  an  rms 
voltage  required  to  attain  90%  of  the  overall  device  transmittance, Vgo-  It  is  unusual  that 
the  turn-off  times  measured  in  our  experiments  vary  as  a  function  of  film  thickness.  The 
turn-off  times  presented  in  Table  9  vary  from  tens  of  milliseconds  for  the  10  fim  shutters 
to  almost  a  full  second  for  the  61  pm  sample.  Vgo  values  presented  in  Table  9  should 
increase  with  increasing  PDLC  film  thickness,  but  instead  are  nearly  constant  for  our 
experiments.  This  indicates  that  the  droplet  size  or  shape  of  the  PDLC  film  varies  with 
thickness. 
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10.2  Transmission  Enhancement 


Measurements  were  also  performed  on  direct  transmittance  of  the  PDLC  films  for  various 
voltages  and  film  characteristics.  These  are  shown  in  Fig.  41. 


Figure  41:  IR  Transmission  as  a  function  of  applied  voltage,  (o)  Germanium  uncoated 
substrate,  d=10  pm  (square)  AR  coated  Ge  substrate,  d=10  pm  (triangle)  AR  coated  Ge 
substrate,  d=26  pm  (*)  AR  coated  Ge  substrate,  d=01  pm. 
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10.3  Angular  Distribution  of  Scattering 


Measurements  of  angular  distribution  of  scattered  light  from  PDLC  films  were  carried 
out. 

Below  is  a  schematic  of  the  experimental  apparatus  used  for  electro-optic  and  differ¬ 
ential  scattering  measurements.  The  in-plane  scattering  plane  is  defined  by  the  incident 
and  scattered  wave  vectors  (k  and  kl  respectively)  and  the  angle  6  between  the  optical 
axis  and  kl  .  The  black  body  radiator  (Infrared  Industries)  was  operated  at  900  de¬ 
grees  C  for  all  measurements.  The  black  body  radiation  was  modulated  at  1  kHs  and 
the  beam  collimated  by  a  zinc-selenide  lens  and  aperture  combination  on  a  layered  InSb 
(2-5  pm)/MCT  (8  -  14  pm)  detector  configuration  contained  in  a  dewax  mount  (E,  G 
k  G  Judson).  The  detected  signal  was  amplified  by  a  dual  preamplifier  and  processed 
by  a  lock-in  amplifier  (Stanford  Research  Systems  Model  SR510).  The  output  of  the 
lock-in  amplifier  was  fed  to  a  Zenith  PC  for  scattering  measurements  and  an  oscilloscope 
(Hewlett-Packard  Model  54501  A)  for  measuring  electro-optic  properties  (transmittance 
in  the  ON  and  OFF  states,  driving  voltage  and  turn-off  time).  Infrared  PDLC  shutters 
were  modulated  with  a  gated  AC  pulse  of  varying  amplitude  and  duration.  FTIR  spectra 
were  recorded  with  a  Perkin- Elmer  Model  1310  IR  spectrophotometer. 

The  following  figure  illustrates  the  effect  of  droplet  sue  on  angular  scattering.  The 
PDLC  was  26pm  thick  with  anti-reflective  coated  Ge  substrates.  The  incident  radiation 
was  polychromatic  and  unpolarised.  As  can  be  seen  from  the  figure,  the  effect  of  droplet 
size  on  the  scattering  profiles  is  pronounced.  For  2-5  pm  and  8  -  14pm  incident  radiation, 
there  is  an  increase  in  the  angular  scattering  with  droplet  size,  as  expected.  Very  little 
scattered  radiation  is  detected  for  8-14  pm  radiation  when  the  droplet  size  is  «  2  pm. 
It  should  be  noted  that  in  both  cases,  only  the  detected  intensity  at  6  >  13  degrees  is 
truly  scattered  and  transmitted  radiation,  evidenced  by  comparing  the  background  curves 
to  those  for  the  PDLC.  At  6  =  0  degrees  (normal  incidence),  only  transmitted  radiation 
is  detected.  Differences  in  the  off-state  transmission  at  normal  incidence  may  be  due  to 
back-scattering.  No  back  scattering  was  evident  between  6  =  90  degrees  and  6  =  150 
degrees.  Angles  beyond  150  degrees  were  inaccessible  due  to  geometric  constraints  of  the 
experimental  set-up. 

These  figures  show  that  there  is  very  strong  forward  scattering.  Beyond  8-10  degrees 
intensity  has  dropped  by  over  two  orders  of  magnitude  in  all  cases. 
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SCA1TEKKD  INTENSITY  SCATTERED  INTKNSIT 


Figure  42:  Angle  dependence  of  scattering  from  a  26  pm  thick  AR  coated  Ge  substrate 
PDLC  with  (a)  2-5  pm  and  (b)  8-14  pm  incident  radiation.  Droplet  sises  are  «  2 pm 
(dash-dot),  10  -  25 pm  (dots),  15  —  30 pm  (solid)  and  background  (dashes). 


11  Summary,  Conclusions  and  Recommendations 

11.1  Summary 

This  investigation  has  concluded  that  there  are  a  number  of  large  optical  nonlinearities  in 
liquid  crystals  which  have  potential  application  in  Optical  Power  Limiting  devices.  They 
cover  the  temporal  range  of  CW  to  approximately  one  nanosecond.  We  are  less  certain  of 
the  applicability  of  liquid  crystals  in  the  picosecond  domain,  but  there  is  reason  to  expect 
that  the  useful  nonlinearities  extend  to  subnanoseconds.  We  are  not,  however,  ruling  out 
the  usefulness  in  the  picosecond  regime. 

In  the  scientific  literature,  both  before  and  during  the  course  of  this  investigation, 
there  were  reports  of  various  large  optical  nonlinearities  in  liquid  crystals  over  this  whole 
temporal  domain.  A  fundamental  element  of  our  investigation  was  the  measurement  of  a 
variety  of  liquid  crystals  to  determine  whether  the  reported  nonlinearities  were  peculiar 
to  a  one  or  a  few  liquid  crystal  species  and  whether  larger  nonlinearities  were  available 
through  a  simple  process  of  searching.  A  second  fundamental  element  of  our  research 
was  to  clarify  the  mechanisms  responsible  for  the  observed  nonlinearities.  Finally,  we  also 
participated  in  analysis  of  potential  applications  and  device  prototype  development. 


11.2  Conclusions 


The  technical  data  presented  in  this  report  and  the  Quarterly  Scientific  and  Technical 
Reports  as  well  as  the  information  obtained  from  other  contractors  during  this  period 
lead  us  to  the  following  conclusions. 


•  In  the  CW  regime,  there  is  certainly  sufficiently  effective  OPL  behavior  of  a  variety 
of  liquid  crystals  for  device  development.  The  active  feedback  prototype  delivered 
under  this  contract  clearly  demonstrated  one  particular  design.  This  design  is  ro¬ 
bust,  capable  of  being  fitted  to  a  variety  of  devices  needed  in  the  field  and  capable 
of  offering  protection  on  a  wide  range  of  time  scales. 

The  passive  OPL  behavior  that  we  investigated  in  this  contract  is  much  less  devel¬ 
oped  than  the  active  form.  There  is  need  for  a  fundamental  alteration  of  the  manner 
in  which  a  passive  OPL  device  is  designed  if  it  is  to  have  the  effectiveness  of  the 
active  device.  There  is  reason  to  believe  that  such  alterations  are  possible,  but  they 
will  at  present  have  to  be  tailored  to  the  specific  environment  and  specifications 
of  field  devices.  If  certain  device  specifications  allow  compromise,  such  as  response 
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time,  temperature  range  or  insertion  loss,  then  variations  such  as  dyed  liquid  crys¬ 
tals  or  temperature  controlled  liquid  crystals  may  be  acceptable  in  passive  devices. 
In  any  case,  since  passive  functionality  necessitates  that  the  active  characteristic 
for  power  limiting  be  inherent  to  the  material,  more  research  on  the  mechanisms  of 
alternative  liquid  crystal  materials,  such  as  mixtures  and  composites,  will  be  needed 
to  move  this  form  of  OPL  action  towards  field  applications. 

•  The  regime  of  nanosecond  OPL  was  one  of  the  other  successes  of  this  investiga¬ 
tion.  Our  group  was  the  first  to  offer  5CB  as  an  effective  OPL  material.  This 
initiated  the  discovery  of  several  materials  with  large  nonlinear  optical  parameters. 
(Other  contractors  discovered  other  liquid  crystal  materials  with  similarly  encour¬ 
aging  properties.)  This  made  it  clear  that  there  was  device  potential  in  this  time 
regime  as  well.  A  major  effort  was  undertaken  to  determine  the  mechanism  for  the 
observed  nonlinearity.  This  effort  was  undertaken  to  avoid  random  searching  for 
nonlinear  enhancement  as  well  as  a  way  to  obtain  a  theoretical  understanding  that 
could  be  used  in  device  modeling. 

Our  investigations  into  mechanism  were  successful,  but  left  several  question  incom¬ 
pletely  resolved.  The  positive  aspect  of  our  results  is  that  we  now  know  that  effective 
OPL  is  accompanied  by  more  than  simply  a  large  nonlinear  refractive  index.  We 
found  strong  nonlinear  absorption  and  the  presence  of  higher-order  nonlinearities. 
The  two  photon  absorption,  which  we  and  others  observed,  seems  to  be  an  im¬ 
portant  mechanism  for  removing  beam  intensity  in  the  OPL  process.  While  this  is 
reasonable  in  and  of  itself,  two  other  questions  are  raised.  First,  why  is  the  nonlinear 
absorption  increased  by  the  presence  of  dyes  (linear  absorption)  in  the  liquid  crystal. 
Second,  where  does  the  energy  go  which  has  been  absorbed.  Simple  conversion  of 
electromagnetic  energy  to  thermal  energy  (by  whatever  relaxation  mechanism)  can 
not  be  effective  in  the  nanosecond  time  scale.  The  energy  blocked  by  the  OPL  is 
enough  to  vaporize  the  sample.  This  is  both  a  positive  and  negative.  It  is  positive 
because  other  competing  materials  which  rely  on  absorption  as  the  OPL  mechanism 
will  suffer  heating  and  thus  have  a  restricted  dynamic  range.  Our  damage  studies 
showed  that  the  liquid  crystal  could  not  be  damaged  (instead  the  sample  cell  or 
surface  alignment  layer  was  damaged).  Such  is  not  the  case  for  some  proposed 
semiconductor  based  OPLs.  The  negative  aspect  of  the  role  of  nonlinear  and  linear 
absorption  is  that  we  do  not  finally  know  what  the  operative  absorption  mechanism 
is  or  how  it  enhances  the  overall  nonlinear  activity  of  the  OPL. 

•  The  effectiveness  of  liquid  crystals  as  optical  power  limiters  in  the  picosecond  regime 
is  still  an  open  question.  It  is  clear  that  new  mechanisms  become  operative  on  this 
time  scale,  and  it  was  one  of  the  successes  of  these  investigation  that  we  discovered 
that  the  nonlinear  activity  of  the  liquid  crystal  continues  to  function  when  the  pulse 
duration  is  reduced  by  a  factor  of  a  thousand.  The  reason  for  the  nonlinear  behavior 
of  the  liquid  crystal  in  this  regime  is  almost  certainly  contains  a  contribution  from 
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purely  electronic  hyperpolarizability  of  the  molecules.  If  this  is  indeed  the  primary 
effect,  then  it  is  likely  that  nonliquid  crystalline  organic  materials  (e.g.  highly 
conjugated  polymers)  will  be  found  to  be  superior  in  OPL  applications. 

Still,  liquid  crystals  offer  potential  in  the  picosecond  regime.  First,  it  is  possible  that 
mechanisms  in  addition  to  the  electronic  mechanism  are  operative.  Such  would  not 
be  expected  in  the  case  of  conjugated  polymers.  The  other  mechanisms,  such  as  long 
axis  rotation,  go  hand  in  hand  with  other  normal  physical  properties  such  as  the 
fluidity  of  the  liquid  crystals.  Fluidity,  for  example,  also  conveys  another  desirable 
effects,  that  of  self  healing  in  the  event  of  damage.  (If  a  OPL  medium  is  fluid,  then 
a  damaged  area  will  disperse.) 

•  The  IR  Shutter  project  was  funded  as  part  of  this  contract  even  though  it  had  noth¬ 
ing  to  do  with  nonlinear  optics.  It  was  however  based  on  liquid  crystal  composites 
and  it.  was  relevant  to  the  mission  of  the  funding  agency. 

While  this  research  project  established  the  shuttering  capability  of  liquid  crystal 
composites  in  the  IR,  it  fell  short  of  establishing  the  materials  as  a  potential  re¬ 
placement  of  mechanical  shutters  presently  in  use.  This  project  did  not  conclude 
that  liquid  crystals  were  unsuitable  for  shuttering  applications,  it  simply  did  not 
establish  that  they  were  suitable. 

When  the  project  ended,  the  use  of  a  liquid  crystal  shutter  in  an  IR  imaging  system 
had  been  demonstrated.  We  were  not  able  to  show  the  potential  use  for  in  night 
vision  applications  and  did  not  find  materials  that  fully  met  the  specifications  we 
were  given  for  night  vision  applications. 

11.3  Recommendations 

In  all  the  areas  described  above,  considerable  experimental  and  theoretical  under¬ 
standing  of  the  systems  was  achieved.  The  value  of  this  data  would  have  been 
greater  had  a  better  understanding  of  the  interface  of  the  contract  work  and  the 
system  applications  been  achieved  earlier  in  the  project.  For  example,  we  were  three 
years  into  the  project  before  written  criteria  for  measurement  of  OPL  performance 
were  provided  to  the  contractors,  and  it  was  a  year  after  that  before  the  first  in¬ 
dependent  measurements  in  an  actual  device  were  performed  by  the  contracting 
agency.  Even  then,  we  did  not  receive  feedback  on  the  performance  of  the  device 
prototype  that  we  provided. 

The  continuation  of  work  in  this  area,  the  use  of  liquid  crystals  in  optical  power 
limiting  applications,  is  strongly  advised.  As  a  result  of  the  past  research,  enormous 
advances  toward  that  application,  as  well  as  scientific  advances,  have  been  made. 
Our  recommendation  is  that  future  efforts  involve  a  more  informed  cooperative 
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venture,  where  the  researcher  does  research  and  the  device  developer  addresses 
engineering  issues  while  communicating  with  each  other. 

There  are  really  two  avenues  that  are  separately  identifiable  as  needed  for  advance¬ 
ments  that  are  not  just  hit  or  miss.  These  are  further  measurement  and  mechanism 
study.  First,  having  identified  the  dozens  of  other  compounds  that  maybe  useful, 
they  should  be  measured  for  the  various  nonlinear  parameters  and  tested  for  optical 
power  limiting.  This  is  more  than  just  putting  it  in  a  sample  holder  and  seeing 
the  power  limiting  curves.  There  is  clearly  more  going  on  than  just  a  bigger  and 
better  nonlinear  index  of  refraction,  so  both  power  limiting  properties  and  nonlinear 
parameters  need  to  be  measured  for  about  30  or  40  other  compounds.  Secondly, 
the  mechanism  for  the  power  limiting  is  still  obscure  and  needs  more  study.  We 
have  managed  to  eliminate  more  wrong  guesses  than  we  have  been  able  to  identify 
mechanisms.  No  device  is  optimised  during  its  development  and  the  devices  under 
consideration  here  can  benefit  greatly  from  a  well  defined  program  that  builds  one 
what  we  learned  already  in  materials  and  device  physics. 
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A  Appendix:  Materials  Survey 


To  help  separate  and  identify  the  operative  mechanisms  for  OPL,  we  systemised  our 
selection  of  the  class  of  materials  to  be  studied.  The  following  properties  and  selection 
criteria  were  applied  in  our  literature  search. 


1.  A  solid  to  nematic  or  smectic  transition  should  occur  near  or  below  room  temper¬ 
ature  and  have  a  range  of  at  least  10  degrees  C  or  more  to  facilitate  experimental 
measurement. 

2.  The  material  should  be  chemically  stable  to  eliminate  uncertainties  about  contami¬ 
nates  contributing  to  or  masking  nonlinear  effects.  Similarly,  laser  induced  reactions 
would  unnecessarily  complicate  analysis  at  this  stage. 

3.  Transparency  in  the  IR  and  visible  range  is  highly  desirable. 

4.  Chemical  breadth  was  desired.  We  sought  compounds  with: 

(a)  Varying  sizes  of  conjugated  cores 

(b)  Varying  linkages  between  benzene  rings 

(c)  Varying  substituent  groups 

(d)  Ease  of  synthesis  or  even  possible  commercial  availability. 

(e)  Simple  and  small  conjugation  structures  to  facilitate  theoretical  calculations 


Except  for  requirements  that  a  nematic  or  smectic  exist  at  a  temperature  accessible  to  the 
laser  measurement  facility,  none  of  the  above  criteria  were  deemed  absolute.  We  compiled 
a  list  of  materials  which  satisfy  all  or  most  of  the  above.  This  list  was  used  by  us  in  the 
selection  of  nonlinear  materials.  Those  selected  were  purchased  or  synthesized  and  are 
listed  in  Sec.  8.2  along  with  additional  materials  used  in  the  study.  This  list  is  provided 
as  part  of  our  recommendations  for  further  materials  evaluation.  Appendix  A  and  the 
tables  in  Sec.  5.1  contain  the  measurement  results. 


108 


AC np  Npo  Is  ->1013 


1. 


I 

I 


CnH2n+  1CO2 


— ^^~^^--CmH2n+ 1 


m  K  -  N  -  I 
6  36  37.5 

6  38.5  44 

.  Chem  15, 441  (1975),  Zaschkett,  H. ,  Stolle,  R. 


CnH2n*l  S 


=  alkyl,  alkoxy  n  K  ♦  8a  ♦  I 

4-6  -25  -76 

.Chem.  17,293(1977). 


CnH2n ♦  1 


O 


CmH2m  ♦  t 


m  K  ♦  N  ♦  I 
4  36  (.32) 

8  35  60 

|  J.  Prakt.  Chem.  317, 617  (1973),  Zaschkett,  H. 

I 

I 

I 


6 


d  <N  co  ■»# 


4. 


CnH2n+l  C6Hl3 


K  -  Sa  -  I 
30  (.28) 

23  32 

24.3  44 


n=8 


K  -  SA  -  N  -  I 
37.4  47.8  86.9 


6. 


C6H9-CH-CH2 

I 

CH3 


CN 

-0 — 0-coo  -h-  Cn\H2m  +  l 


m  K  -  N  -  I 

8  29  71 

5  37  49.5 

CA  88, 137448 

Erdmann  D.,  DZ-OS  2613293  (1977) 


7 


7. 


CnH2n+  1*C*C 
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n  K  ♦  N  ♦  I 

7  23  79 

8  27  81.5 

9  36  80 

Angew.  Chem.  86, 378  (1976). 
Van  Der  Veen,J. 


C„H2n+i  -(Q^-  CH=N-N=CH-^^-CmH2m+i 


n  m 

K  ♦  N  ♦  I 

-  ^  - T-  4 

3  6 

23  92 

3  8 

28  89 

Mailer,  H.,  DE-OS  2852881  (1979). 

• 

MCLC  35, 155  (1976),  MCLC  41, 141  (1978). 
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13. 


C2H5O  — CH = N  — Cg2 •  CH- C2H5 

CH3 


K  *  N  -  I 
29  60.1 

J.  Chtm.  Phys.  58, 413  (1973) 


14. 


CnH2n+l 


CH=N 


CmH2m  +  l 


n  m  K  *•  Sb  ♦  I 
3  5  27.1  37.1 


15. 


C6Hi3 

C6H13 


K  ♦  Sb  N  ♦  I 
20  44 

RaCN  K  ♦  N  ♦  I  CN 

30 

MCLC  67, 241(1981). 
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16. 


Cl 


K  ♦  N  ♦  I 
35 

ABS  4,  LCCSL 1981  (Cll) 


n  m  K  ♦  Sa  *  N  ♦  I 

8  4  22.5  31.5  33 

International  Liquid  Crystal  Conference ,397 ,  (1973) 
Am.  Phys.  (Paris)  3. 131  (1978). 


CnH2n+l— ^D)-CH=N  — 


n  K  N  I 
6  26 

7  25  (8.5) 

Billard,  J.,  Dubois,.  J.C.,  J.  Phys.  (Paris)  Coll.  Cl  36,  Cl-355  (1975). 
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19. 


Cm  H2m+1 

— ^^~C02 

— O2C  — CmH2m+l 

Cn  H2n+I 

m  n 

K  ^ 

N  ♦  I 

7  6 

19 

41 

5  6 

6  6 

33.5 

29 

41 

Crys.  Res.  Technol.  19, 55  (1984),  Weissflog,  W„  Demus,  D. 


20. 

CnH2n+l  — ^3V~COS 

n  m  K  N  ♦  I 

5  4  28  34.7 

4  4  13.9  22.5 

Reynolds,  R.,  et  el.,  MCLC  36, 41  (1976). 


21. 


CnH2n+l  CO2  — 0"CmH2ati 


T 


n  m  K  -  N  -  I 

Several  -room  temp.  -R.T. 

MCLC  67  1(1981). 
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22. 


C02— ^-R 

N  -  I 

.2.0  r.t. 

CO 

50 

MCLC  67, 1  (1981). 


a  R 

1  CH3 

to  F 


23. 


C2  H5  -  CH  -  CH2  C02  CN 

ch3 


K  ♦  N*  -  I 
37 

J.ACS  90, 1585  (1985). 


24. 

c5  Hu  H^)-  co2  -<^-CH0 


K  *  N  -  I 
33.8  (.24) 

Osman,  M.A.,  DE-OS2716610  (1978). 


13 


25. 


Cm  CO2  — CnH2n+l 


m 

1 


42 


Z  Naturtosch  27b.  774  ( 1972),  Steinstraser,  R. 
Arora,  S.,  U.S.  -  Pat  4086002  (1978). 


26. 

Cm  H2m  +  1  co2 


n  m 

vary 


K  -  N  ♦  I 
15Q-50#  15“-50* 


27. 


C2H5-CH-CH2 

CH3 


OC„H2n. 


n  K  *  N  ♦  I 

5  26  42 

6  30.5  51.8 

CA  95, 124545  (1981) 


14 


n  m  K  *  N  -  I 

6  5  29  48.3 

Zh.  ObsLch.  Khim.  46, 2125(1976),  Daugvila,  J. 


29. 

Cn  H2n  +  1  *X 


Cm  H2m-i 


n  m  K  -  N  -  I 

50-60  X=_v  = 

vary  RT 

MCLC  75, 95  (1981);  67, 109  (1981). 


30. 


CnH2n-rl 


CN 


n  K  -  N  -  I 

4  31  44 

CA  92, 6287(1978) 


-.O.CO 


15 


31. 


C„Hj„4.|  -OcOS  m  Him  ♦  l 


K  -  N  -  I 
—20  37-40 


Sato,  H.,  et  al.,  CA  91. 210978  (1979). 
32. 


22.6  38.7 

Helv.  Chim.  Acta  65, 1242(1982),  Petrzilka,  M. 


33. 


CnH2n-H— (  C02  — ^Q)-CmH2m+l 


K  N  I 

R.T.  R.T.or— 50° 

Gray,  G.,  et  al.,  MCLC  67, 1  (1981). 

CnH2n*i— ^  )—  00,  )— CmH2m-t 

K  -  N  -  I 
R.T.  35-3CT 

Adv.  in  LC  Research  (1980),  S1075,  Deutscher,  H.  J. 
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A  Appendix:  Data  Sheets  and  Results  Tabulation 


The  attached  data  sheets  and  tabulation  of  properties,  including  x3  values  and  nonlinear 
parameters  were  provided  as  part  of  the  Scientific  and  Technical  Reports  during  the  course 
of  the  contract  and  are  reproduced  here  for  completeness. 


119 


tame  of  OipmittlKM:  LIQUID  CRYSTAL  INSTITUTE,  KENT  STATE  UNIVERSITY 

Puncyai  tawaigawr  Michael  A.  Laa  (Co-Pie:  Peter  Palf  fy-Muhorny ,  John  L.  Weet,  Liang-Chy  Chian,  J.  Will  tan  Soane! 


«  8S*C  for  isotropic  phase 


Name  of  LIQUID  CRYSTAL  INSTITUTE.  KENT  MAIL  UHl  vut..i  •  • 

Principal  Investigator:  Michael  A.  Lee  (Co-PIs:  Peter  Pnlf fy-Muhorny ,  John  L.  West.  Llang-Cliy  Chien.  J.  Wllliaw  Doane) 


T=97*C  for  isotropic  phase 


MAlliKIALS  PROPERTIUS  QUESTIONNAIRE 
(See  MstrvciK*  sftcci) 


m 


IKO*C  for  nematic  phase 


(Sec  MStmetMM  shod) 

Name  ol  Oipnimio*.  LIQUID  CRYSTAL  INSTITUTE,  KENT  STATE  UNIVERSITY 

tancioai  Inmlijraoc  Michael  A.  Lee  (Co-PIs:  Peter  Palffy-Muhoray.  John  L.  West,  Liang-Chy  Chien,  J.  Wlllla*  Doana) 


IR0*C  Tor  nematic  phase 


(Sec  MtinClMM  kitcclj 

tiamc  of  Oigmuiioa:  LIQUID  CRYSTAL  INSTITUTE,  KENT  STATE  UNIVERSITY 

lYincipat  famiifMor  Michael  A.  Lee  (Co-PIa:  Peter  Palffy-Muhorey,  John  L.  West,  Liang-Chy  Chian,  J.  Wl  Ilian  Doaa*) 


*N«  measurable  nonlinearity  wu  observed  within  our  experimental  error! 


not  measured  because  our  temperature  cannot  go  below  22 *C 


{See  mmtikik*  vlicci) 

Nwc  of  OtpniiMk*:  LIQUID  CKYSTAL  INSTITUTE,  KENT  STATE  UNIVERSITY 


Arc  for  isotropic  phase  0X  =  30  (£* ) 


MATLKIALSi  PKOPliNTUiS  UWJil  H/NHAIMli 
(Sec  Miiraclio*  ttaci) 

ftawc  trf  Orptimnwr  LIQUID  CRYSTAL  INSTITUTE,  KENT  STATE  UNIVERSITY 

Phac*«l  iHmlifMT  Michael  A.  Lee  (Co-PIs:  Peter  Palf fy-Mul.oray ,  John  L.  West,  Liang-Chy  Chian.  J.  Witllan  boane) 


80*C  for  isotropic  phase  p.  _  n(^) 


(.S'cc  HIMMKIWH  tiKC I) 

Nm,,c  of  OigjniiMHm:  LIQUID  CRYSTAL  INSTITUTE,  KENT  STATE  UNIVERSITY 

Fhncipil  iHvesliptor  Michael  A.  Lee  (Co-PIss  Peter  Palffy-Muhoray .  John  L.  West,  Ltaug-Chy  Chien,  J-  Wliliaei  Doane) 


.  Name  of  Oigaiizalioa:  LIQUID  CRYSTAL  INSTITUTE,  KENT  STATE  UNIVERSITY 
fYincipal  l«.veai|*oc  Michael  A.  Lee  (Co-PIs:  Peter  Palf  fy-Muhoray .  John  L.  West,  Liang-Chy  Chien.  J.  WllUa.  Doane) 


not  measured 


